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Abstract 



This thesis is organized in two independent parts, in which I study two different aspects of high 
T c superconductivity. 

The first part begins with an introduction aimed to briefly introduce some relevant exper- 
imental and theoretical works performed in recent years, that have helped us to think about 
cuprates the way we do now. Afterwards, I introduce Landau Fermi liquid theory in a stan- 
dard text book way. The question of validity of Fermi liquid theory in 2-d is then raised and 
investigated by searching for singularities in Landau's /-function. I show that the interaction 
function between two quasiparticles whose momenta approach each other near a curved point 
of the Fermi surface, contains a 1-d singularity not strong enough to change the Fermi liquid 
behavior. On the other hand, inflection points provide 2-d singularities that have to be taken 
seriously in Fermi liquid considerations. I then introduce nearly antiferromagnetic Fermi liquid 
theory (NAFL), which is a phenomenological theory proposed to describe high T c systems. I 
mainly focus on the self-consistency of the theory in calculations. I criticize the theory on the 
basis of overlooking the vertex corrections in the strong coupling calculations of the transition 
temperature T c . I calculate the first vertex correction for an optimally doped system and show 
that it is of the same order of magnitude as the bare vertex. Migdal's theorem is therefore 
not valid and Eliashberg formalism is not applicable to this situation. The same conclusion is 
obtained even after inclusion of the quasiparticle residue Z to the calculation. The sign of the 
vertex correction is then considered. I show that the positive sign of the vertex correction for 
the optimally doped system requires a phase transition of some sort as the doping is decreased. 

Part II of the thesis is devoted to the vortex lattice properties of high T c superconductors. 
I establish a method to study vortex lattice properties of d-wave superconductors based on a 
generalization of the London model. The method has the advantage of simplicity as well as 



ii 



having very few free parameters (one at most) compared to other methods. The generalized 
London free energy is obtained from an s-d mixing Ginsburg-Landau free energy and also from 
the microscopic theory of Gorkov. The generalized London equation is found to be analytic 
at high temperatures. At very low temperatures however non-analyticities arise as a result of 
the nodes on the superconducting gap. I then present the results of our calculations of some 
measurable quantities, such as the vortex lattice geometry and the effective penetration depth 
(as defined in /j,SR experiments). Comparison between our results and different experimental 
data is then performed. Especially our prediction for the magnetic field dependence of the 
effective penetration depth at T = 0, which is recently observed in /xSR experiments with 
excellent agreement, is discussed. 
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Chapter 1 



Introduction 



It has been more than a decade since the discovery of the first high temperature superconductor 
0], and still no consensus about the microscopic theory relevant to cuprates exists. New theories 
have been emerging one after another and disputes about the older theories are still unsettled. 
A general feature of all cuprates is their strong anisotropy which makes them effectively 2-d 



electronic systems. Fig. shows the crystal structure of the YBCO compound which is one of 
important high T c materials. At the center there are two Cu02 planes separated with Y atoms 
in between. Between these two planes and another two planes above or below them, there 
are oxygen chains (that play crucial role in the doping of the system) and also Ba atoms. It is 
generally agreed that the Cu02 planes play an important role in determining the characteristics 
of cuprates at all dopings. Fig. [O] shows the crystal structure of YBCO material In the early 
days of high Tc's, Anderson [g] proposed an effective Hamiltonian which describes a single band 
2D Hubbard model with an infinite on site repulsion U. Soon after, Zhang and Rice || tried to 
justify the Anderson's suggested model (known as the t-J model) for cuprates. They proposed 
that the holes which primarily reside on oxygen sites would delocalize onto four O sites around 
Cu +2 ions. The combination of a delocalized O hole and a Cu +2 ion forms a singlet which 
would hop between Cu +2 sites in the same way as a hole does in a simple Cu +2 square lattice. 

In the absence of doping, Cu02 planes are half filled and form a quasi-2D antiferromagnet 
with an insulating Neel ground state. Numerical simulations, experiments, and scaling theory, 
all suggest that above the Neel temperature, the antiferromagnetic spin correlations are well 
described by the isotropic 2D Heisenberg model H with exchange coupling J ~ 1550-fsT between 
nearest neighbors. However, according to the Hohenberg-Mermin- Wagner theorem M, classical 
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Figure 1.1: Crystal structure of YBCO 

fluctuations destroy Neel ordering at an arbitrarily small nonzero temperature. The small 
(~ 10~ 5 J), but finite correlations between the planes are believed to be responsible for the 
observed long-range order [Q]. Chakravarty, Halperin, and Nelson |J on the other hand, have 
employed quantum non- linear a model for the long- wavelength action of the 2D, s = 1/2 
Heisenberg model. A scaling analysis of this model describes correctly both experimental and 
numerical data [Q]. 

Introducing a small amount of holes into the Cu02 planes, the Neel temperature falls off 
rapidly from its undoped value. Increasing the doping level even further, the system moves into 
a state with a metallic normal state and a transition to a superconducting ground state at low 



temperatures (cf. Fig. 1.2). The superconducting transition temperature increases initially as 
doping is increased. After reaching its maximum at optimal doping, it starts to decrease and 
the superconductor enters a so called overdoped regime (Fig. 1.2). All over the metallic phase, 
the antiferromagnetic correlation is short ranged with £ between one and (up to) fifteen lattice 
spacings, depending on the material, doping and temperature. 
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Figure 1.2: Phase diagram of high temperature superconductors, x represents doping, AF 
antiferromagnetically ordered phase and SC superconducting phase. Upper line is the crossover 
to pseudogap state. [^] 

Experiments evidently indicate that the electronic properties of the underdoped and over- 
doped materials are different (7|, ||, ||, 10, O]; especially the pseudogap phase which exits only for 




underdoped compounds |12j (Fig. |l.2|) . Among experiments, inelastic neutron scattering (INS) 
and nuclear magnetic resonance (NMR) experiments provide important information about the 
magnetic properties of the system. INS directly probes the spin response function, w). In 
Lai.85Sro.i5Cu04 the INS experiments |§ display four sharp magnetic peaks located at Sq ~ ir/4 
away from the antiferromagnetic wave vector. These incommensurate peaks were also observed 
recently in YBa2Cu306.6- [0] Like many experiments, INS has some inherent problems. For 
example it is difficult to determine experimentally the absolute value of intensity. It is also not 
easy to separate the phonon spectrum from the magnetic excitations. 



NMR experiments |Tj, 15] on the other hand, are much more exact, and can detect much 
lower intensities. However, NMR is an indirect probe and its data analysis involves using 
a theoretical Hamiltonian connecting nuclear and electron spins, the hyperfine Hamiltonian. 
Moreover, since NMR experiments are done at very low frequencies (~ nuclear spin Zeeman 
splitting), no information about w-dependence of x( w ) T) can be obtained. Early measurements 
of the Magnetic Hyperfine Shift (Knight shift) at the oxygen sites showed that most of the total 
spin polarization was carried by the copper sites [16, O]. Alloul et al. |H| reported that in 
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YBa2Cu306+x) the knight shift of Yttrium scales linearly with the macroscopic susceptibility. 
This was interpreted as a proof of the existence of a single spin fluid in agreement with the t-J 
model, because yttrium is coupled to the magnetic susceptibility of the Cu02 planes mainly 
through the polarization of the oxygen orbitals. Further measurements of knight shift showed 



the same T-dependence on all nuclear sites [14]. Similar result was obtained for all other 
cuprates indicating that the single spin fluid nature is a common feature of all these materials. 
Assuming a single fluid picture, Shastry and also Mila and Rice pOl (SMR) proposed a 
hyperfine Hamiltonian in which the only source of coupling is the localized electronic spins at 
the copper sites. 

While the magnetic hyperfine shift (Knight shift) gives a measure of the uniform magnetic 
susceptibility, the Nuclear Spin Lattice Relaxation Rate (1/Ti) provides information about 
the imaginary part of the dynamic spin susceptibility %(<f,a;). An important experiment by 
Takigawa |l?J] on YBa2Cu307, in the early days, displayed a significant difference between the 
temperature dependence of 1/T\ for the Cu02 plane copper and oxygen sites. The relaxation 
rate for oxygen sites (and also 89 Y sites) obeys Korringa relation 17 TiT ~ Const, while the 
relaxation rate for copper sites is much shorter than for 17 and 89 Y sites and shows non- 
Korringa behavior. Using SMR Hamiltonian, Moriya (21]] showed that 

mr^ochrnVF^^l (1.1) 
q 

where yf'{q,ijS) is the imaginary part of the magnetic susceptibility and F(q) is a form factor 
which depends on the site and the direction of the magnetic field. For 17 sites, F(q) is peaked 
at q = and vanishes at the antiferromagnetic (AF) wave vector Q = (7r/a,7r/a). For 63 Cu sites 
on the other hand, it is peaked at q = Q. Thus different nuclei probe different regions of the 
momentum space of \"{q,ti)- The significant difference between the 63 Cu and 17 relaxation 
rates could then be understood as a result of an enhanced susceptibility at q = Q. 

Another important quantity measured by NMR experiments is the nuclear spin-spin re- 
laxation rate (I/T2G) which is the rate of relaxation of the nuclear spins among themselves 
and is usually much larger than \jT\. In most solids Tig is temperature independent and is 



Chapter 1. Introduction 



5 



related to spin-spin dipolar coupling and therefore doesn't give any information about electron 
susceptibility. In cuprates however, because of the antiferromagnetic enhancement of the spin 
susceptibility, indirect coupling of spins via nonlocal spin susceptibility dominates the nuclear 
spin-spin coupling. Thus T2G can provide important information about the antiferromagnetic 
exchange between electron spins as was pointed out first by Pennington et al. [j22[|. Measure- 
ments of 1 jT<iG in different materials with different dopings have found 1 )Tig to increase with 



decreasing temperature, with sometimes a flattening or small decrease close to T c [14|. This is 
considered as a proof for the existence of AF correlations in all cuprates. 

To have a quantitative description, Millis, Monien and Pines (MMP) (2^| proposed a phe- 
nomenological form for the magnetic susceptibility x(<f) w ) which has a peak near the Antifer- 
romagnetic (AF) wave vector Q = (ir/a,ir/a) 

x(q, w) = ^ + — ^V- (1.2) 

l + i 2 {q-Q) 2 -iuj/usf 1-w/Io 

where xo is the measured bulk spin susceptibility which in general is temperature dependent, 

Xq is the staggered susceptibility which is some orders of magnitude greater than the bulk 

susceptibility and Tq is an energy scale of order of quasiparticle bandwidth or Fermi energy. 

The first term in Eq.(|1.2j) produces a peak at AF wave vector Q and the second term is the 

ordinary Fermi liquid susceptibility. For the 17 sites, the effect of the first term is suppressed 



by the form factor F(q) (cf. Eq. (LI)) and the second term (the Fermi liquid term) is the only 
term that remains. To the lowest order, Eq.( |1.2| ) gives ( 17 TiT) -1 oc xo(^) which is consistent 
with the Korringa relation. For copper sites however, the contribution of the first term is 
dominant near Q and one can safely neglect the second term. Assuming large correlation 
length £, the form factor can be treated as a constant and the integration can be simply carried 
out to get 

( 63 TiT)~ l oc -?%r = — (1.3) 

where a is defined by xq = a £ 2 an d assumed to be temperature independent. Thus the spin 
lattice relaxation rate provides information about the temperature dependence of u s f in the 
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MMP model 

63 IiT oc u sf /a (1.4) 

This is completely different from the 17 case and has non-Korringa behavior. 

In the early days of high T c 's, there existed some speculations that antiferromagnetic cor- 



relation was responsible for Cooper pairing in high T c materials [24, |25|, ^7|, 28| (the idea 
of d-wave pairing already existed in the context of heavy fermions |2^, |30|). Inspired by para- 
magnon theory, Pines and collaborators |HJ, |3^[ proposed a (phenomenological) nearly anti- 
ferromagnetic Fermi liquid (NAFL) theory to describe pairing process as well as normal state 
properties of high T c superconductors. The theory assumes that the low lying excitations of a 
high T c system in its normal state, are Fermi liquid quasiparticles interacting with each other 



via antiferromagnetic spin- fluctuations. The susceptibility (1.2) can be thought of as the prop- 
agator of bosonic particles that describe spin fluctuations. One of the first predictions of NAFL 



theory was d x 2_„2-wave pairing which was then confirmed by various experiments |33|, 34]. 



Transition temperature was calculated in weak [31] and strong [35, [3(| coupling limits and a T c 
as high as 90K was obtained using several fitting parameters. 

NAFL theory, as we mentioned above, is constructed based on the assumption that the 
excitations of the system are Fermi liquid quasiparticles. Almost none of the normal state 
properties of cuprates agree with Fermi liquid theory (at least in low to moderate doping 
regimes). The proximity to Mott-insulating antiferromagnetic phase might suggest that non- 
Fermi liquid properties are actually signatures of non-Fermi liquid excitations. NAFL theory 
on the other hand believes that the anomalous behavior is a result of strong interaction between 
the Fermi liquid quasiparticles due to antiferromagnetic spin fluctuations. 



Orthogonal to the above philosophy was Anderson's proposal in 1990 [37, Q that 2- 
dimensional Fermi liquids are unstable to the formation of a 2-d Luttinger liquid. Since Ander- 
son, theorists have been trying to understand if and when "non-Fermi liquid" (NFL) can appear 



in itinerant Fermi systems. There is strong evidence [39] for NFL behavior in the v = 1/2 quan- 
tum Hall liquid; theoretical understanding of this assumes a singular gauge interaction between 
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the quasiparticles ppj ]. The high T c problem is more complex because the lattice introduces an 
on-site constraint for strong Hubbard repulsion- which means the lattice problem is fundamen- 
tally different from a continuum model. Perturbative investigations of the nature and existence 



of singular interactions in the low-density regime [41, 42, [43|] have not shown NFL behavior, 
even in a finite box |4*|] (although they do reveal weak singular structure in the quasiparticle 
interactions [^, |4^]). However, these results do not rule out non-perturbative singular behavior 
in the higher density lattice problem- no consensus yet exists on whether these exist in 2 dimen- 
sions. We stress that this question is quite different from that of the existence of instabilities 
to ordered phases (e.g., Superconducting or SDW instabilities), which can be understood using 
standard methods Q|. 

As discussed by various authors 37, |3^, NFL behavior is typically expected if the 
interaction between quasiparticles with momenta p and p 1 is singular when Q = p — p' — > 0. 
It is nevertheless essential to satisfy all relevant Ward identities and crossing symmetry in any 
calculation, given the way singular behavior in one channel can influence that in the others. 
An approach which reveals the structure in the irreducible interactions as a function of Q 
and which satisfies the above requirements, involves calculating the irreducible 4-point vertex 
fpp,' perturbatively JHI , and then using this to determine the full scattering function and 
quasiparticle dispersion. If f££,' is not singular, then it is nothing but Landau's /-function. 

In the next chapter, after an overview on Landau Fermi liquid theory, we perform a calcula- 
tion of fpp,' along the line mentioned above for a regular Fermi surface and also near inflection 
points of the Fermi surface. Our calculation for a curved Fermi surface reproduces the previous 
results |43| existing in the literature. Near inflection points however, we find that f££ is singu- 
lar in Q already at 2nd order in a short-ranged repulsive interaction. We stress that the NFL 
behavior found previously for Fermi surfaces having finite curvature everywhere appeared only 
at infinite order, and was too weak to destroy Fermi liquid behavior [[4^, ^3jj . Calculating the 
imaginary part of the self energy to second order in the singular interaction, we find a frequency 
dependence different from the expected Fermi liquid form. 
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Besides the controversy about the Fermi liquid starting point, the validity of Eliashberg 
formalism used in NAFL theory to calculate T c is under question. Migdal's theorem is an 
essential ingredient of Eliashberg formalism which states that the vertex corrections to the 
bare coupling between electrons and the bosonic medium (here spin fluctuations) are small. 
In the case of phonons, the small ratio of Debye frequency and Fermi energy (lud/E f ) makes 



the vertex corrections small [46]; despite the large bare coupling. In the case of spin-fermion 
coupling on the other hand, the only relevant energy scales are E F and the antiferromagnetic 
exchange coupling J. The ratio J/E F however, is of order one |47|| . Thus it is not clear why 
vertex corrections should be small in the spin fermion interaction. 



Calculations of the vertex corrections close to the antiferromagnetic ordered state [47, [4^, [49[ 
showed that the vertex corrections are actually the same order as the bare vertex and therefore 



are not small. Schrieffer [47] showed that in the ordered state, the dressed vertex should vanish 
at the antiferromagnetic wave vector. The vertex corrections therefore should be equal to the 
bare vertex but with opposite sign. The magnitude of the vertex correction at higher dopings 
is more important because the superconducting transition is suppressed at low doping. 

In chapter 3, we first give a brief introduction to paramagnon theory, which is a theory for 
Fermi liquid systems near a ferromagnetic instability. We then turn to the nearly antiferro- 
magnetic Fermi liquid theory pointing out the similarities and differences between NAFL and 
paramagnon theory. The self-consistency of NAFL theory is then discussed. We calculate the 
first correction to the bare spin-fermion vertex. We find that the vertex correction is actually 
big and of the same order as the bare vertex. The self-consistency of the calculation including 
quasiparticle renormalization is discussed afterwards. 

The sign of the vertex corrections is another subject of dispute. As we mentioned above, 
the sign of the vertex corrections is negative near the ordered state. At optimal doping the sign 
of our vertex correction is also negative in disagreement with some other results |5(J. We talk 
about the sign of the vertex corrections at the end of the 3rd chapter. Our discussion is based 
on the flow of the renormalized vertex versus change in doping. We argue that the positive sign 
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of the vertex corrections at the optimal doping requires a phase transition at some intermediate 
doping. In fact, for other reasons, this phase transition has been suggested by Chubukov to be a 



Lifshitz transition in the topology of the Fermi surface [51]. However, there is no experimental 
evidence for this or any other phase transitions at intermediate doping at the time of writing 



Chapter 2 



Fermi Liquid Theory 



2.1 Landau Theory of Fermi Liquids 

Systems of strongly interacting electrons are very difficult, if not impossible, to solve. Perturba- 
tive methods fail to work as a result of the strong coupling between particles. Non-perturbative 
methods, also, are not quite well developed for systems of more than one dimension. It was 
Landau's ingenuity that realized the low energy excitations of a strongly interacting fermion 
system could still be fermionic quasiparticles similar to non-interacting electrons. Proposed 
in 1957 [^, |54j, Landau's Fermi liquid theory has been extremely successful in describing low 
temperature properties of liquid 3 He [55, 56] as well as interacting electrons in metals [57]. It 
especially explains the success of non-interacting electron theory for low temperature properties 
of metals. 

In his theory, Landau assumes that if one adiabatically turns on the interaction between 
electrons, the non-interacting ground state evolves into an interacting ground state with a one 
to one correspondence between the bare particle states and the quasiparticles in the interacting 
state. A necessary condition for this to happen is that no bound states should form during this 
process. A system that possesses this property is called "Normal Fermi Liquid" . Superconduc- 
tors therefore do not belong to this category. In Landau's picture, the number of quasiparticles, 
N, is the same as the number of bare electrons. They carry the same charge as electrons and 
obey Fermi statistics. One has to be careful not to go too far with this picture. For example, 
the ground state energy is not the same as the sum of the quasiparticles energies. As we will see 
later, the quasiparticle states are meaningful only near the Fermi surface. A better way, is to 
define the quasiparticles as excitations of the system above its ground state and also quasiholes 
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as absence of particles in the ground state. More precisely, a quasiparticle state with momen- 
tum p (> p F , Fermi momentum) is defined by adiabatic evolution of a non-interacting Fermi 
sea with an extra electron at momentum p (and likewise a quasihole state). 

Yet another assumption is necessary to establish the Landau Fermi liquid theory. The 
assumption is that the quasiparticles interact with each other via an interaction function fZ~, • 
(For simplicity from now on we suppress the spin indices and denote the /-function by f pp i . We 
put them back when necessary later.) The scattering of quasiparticles due to this interaction 
causes the quasiparticles to decay. If the decay rate is large, then the quasiparticle lifetime, 
t p may be too short. The quasiparticle states therefore will be ill-defined. In other words, the 
quasiparticle states are not exact eigenstates of the interacting Hamiltonian, i.e. the energy 
levels have finite width (~ 1/t p ). In order to have distinguished energy levels, the width should 
be smaller than the energies, or 

- « (e p - E F ) (2.5) 

Tp 

The key point here is that the scattering of these quasiparticles occurs at a vanishingly small 
rate as the quasiparticles get close to the Fermi surface. The reason behind this lies in the Pauli 
exclusion principle which restricts the phase space available for quasiparticles to scatter to. At 
T = this phase space is proportional to (p — p F ) 2 where p is the momentum of the scattered 
particle ]58|] . The scattering therefore, is very small when p ~ p F . At higher temperatures, the 
scattering rate increases like T 2 . As a result, the quasiparticle picture is meaningful only at low 
temperatures and for quasiparticles with momenta close to the Fermi surface. Thus Landau's 
picture, with N quasiparticles filling up the states up to the Fermi surface, should be considered 
only as a formal assumption. Still this picture is useful especially in the context of Luttinger's 
theorem that allows to state that for a spherical Fermi surface, the value of Kp is the same in 
the interacting and non-interacting system. 

Let us apply a weak perturbation to the system to take it away from its ground state. The 
effect of this perturbation is to change the occupation number by 5n^. Landau postulated that 
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the change in the energy of the system is given by the expansion 



1 



SE = ]T epn'l + — ]T frtSrfynl 



(2.6) 



p,p 



where eji is the quasiparticle energy and V is the volume of the system. In general, one can 
expand to higher powers of 5n^. However, expansion to the second order is enough for most 
practical purposes. The quasiparticle energy affected by other quasiparticles is given by 



6 P ~ £ n — e p + y JpP' 0n pi 
p w 



(2.7) 



Landau's /-function is therefore the second derivative of the energy with respect to 8np and 

(2.8) 



Sn -, 
p' 



fpp '~ 5n%5n% 
p p' 



Since there is a one to one correspondence between quasiparticle states and non-interacting 
electron states, the quasiparticles should obey Fermi statistics (a more rigorous proof involves 
maximization of the entropy. (57|) 

1 



(2.J 



p e P{tp-v) _)_ i 

where \x is the chemical potential; which can be shown to be equal to the Fermi energy Ep at 
T=0 |59[| . For a SU(2) symmetric system, one can write (putting back spin indices again) 



pp' J pp ^pp' 



These functions are usually expressed in terms of spherical harmonics 

'21 + V 



(2.10) 



fl 



Ait 
21+1 
Att 



dnp^cosd) fa 



dflPi(cos i 



pp 



\p\=\pi\=p F 



\P\=\p'\=pi 



(2.11) 



fl and (pi are known as Landau parameters. In terms of these parameters, the quasiparticle's 
effective mass is given by (for all derivations and proofs see Ref. f60f) 
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Observable quantities are usually given in terms of dimensionless Landau parameters F[ = 
N(0)fi and Z\ = N(0)<pi, with the density of states at the Fermi surface 

Ar (°) = ^E^-^) = ! ^ (2-13) 



In terms of i*)'s and Zi's, we have 



m* ( F-\ 
Effective Mass : — = 1 + — 



m \ 3 
Specific Heat : cy = -m*p F k 2 B T 

o 

o 

I Tip 

Compressibility : — = — — (1 + Fq) 

k 3m* 

v 2 

Sound Velocity : c? = — ^(1 + F 

3mm* 

2 * 

Magnetic Susceptibility : Xa ~ 



4vr2(l + Z ) 

where n = N/V is the density of electrons and 7 = eh/mc is the electron's gyromagnetic ratio. 
2.2 Microscopic Foundation of Fermi Liquid Theory 

As emphasized in the previous section, the basic ingredient of Landau's Fermi liquid theory 
is the function f pp t; which describes the interaction of low energy quasiparticles. The simple 
form of the interaction term in ( |2,6[) , essential in Fermi liquid calculations, is far from obvious. 
Renormalization group theory jl4| provides a plausible derivation for a low energy effective 
Hamiltonian in agreement with Fermi liquid theory. As we will see later however, renormaliza- 
tion group might miss some essential singularities in the quasiparticle interaction by putting all 
the momenta on the Fermi surface. Moreover, although renormalization group is a conserving 
approximation pl|, B^], it does not satisfy the Pauli principle [63, 34. 65[|. Another way, is to find 



a relation between f pp i and the dressed two-particle interaction or 4-point vertex function in 
many-particle theory. What we do in this section is to first derive f pp i for a low density electron 
gas by variation of the ground state energy using second order perturbation theory. We then 
discuss the connection between many-particle theory and Fermi liquid theory. In particular, 
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we find a relation between Landau's /-function and 4-point vertex function in many-particle 
theory. 



2.2.1 Calculation of f pp / for a Dilute Fermi Gas 



In a dilute Fermi gas, one can calculate f pp i using Abrikosov-Khalatnikov's formalism |45| (for 



more detailed derivation, see also Ref. |66|], or Ref. [43, 67] for derivation in 2-d). In the dilute 
limit one can assume short range interaction between electrons. The interaction Hamiltonian 
therefore can be approximated by 

#int = y E 4+^-^1^ (2 - 14) 

q,p,p' 

where U is the coupling constant and V is the volume of the system. We do our expansion in 
powers of the s-wave scattering amplitude, a, which is given to the first order in U by 

4-7ra 

U = (2.15) 

m 

Applying perturbation theory to H{ nt , the first order correction to the ground state energy is 
given by 

£« = ^5>n 2 Q 12 = ^ (2.16) 
V f£ mv 

Here the index (1) represents both momentum and spin, m is the occupation number and Qij 
requires spins to be antiparallel 

Q %3 = \{\-B l -a j ) (2.17) 

where cr's are the Pauli matrices. To second order in perturbation theory we have 

p( 2 ) l(^int)mo| 2 ^ , niw 2 (l ~ n 3 )(l - n 4 )Qi 2 Q 3 4 ,„ 



Eq. ( 2.1S| ) is divergent at large momenta. To get away from this divergence we have to expand 



everything in terms of a instead of U. First notice that ( |2.15| ) is not exact. To second order in 
U, ( [2.15D can be written as 

^ = U = U + ^T5 + + 9* (2.19) 



Chapter 2. Fermi Liquid Theory 



15 



Here, U is the renormalized coupling. Inverting (2.19), U can be written in terms of a 

ina (47ra) 2 ^ r Q 3 



U 



rn 



y $ 



m 2y Z-^i p 1 +p 2 'P 3 +pa 
1,2,3,4 



(2.20) 



e 1 + e 2 - e 3 - e 4 

substituting this into (2.16), we get a term second order in a which should be combined with 



(|2.1Sj ) to get the second order correction to the ground state energy. Combining all, after some 
manipulations, we get 



2^2 



32ir z a 



m: 



1,2,3,4 



p 2 .P3+P4 



nin 2 ra 3 Qi 2 Q34 
ei + e 2 - e 3 - e 4 



(2.21) 



This now converges at large momenta, since all momenta are limited because of the n^'s. Taking 
the derivatives with respect to n pa and n p ' a >, we get 



PP' 



Ana 



rn 



-$cr,-cr> 



(47ra) s 
m 2 V 



E 



25 n 



^p' ^p+p'—k 



+ - 



1 



+ 



1 



(2.22) 



e P V' + € k £ k+p-p' V' 6 P € k e k+p'-p 

The first term in the bracket is the "Cooper" channel and the last two terms belong to the 
"crossed" channels. The denominators of the last two terms vanish when p = p' . The singu- 
larities however get canceled between the two terms except for special cases (we will return to 
this point later). 



2.2.2 Fermi Liquid Theory and Many- Particle Theory 



In many-particle theory [58, 60, 66, 65], the excitations of the system above its ground state 
are given by the poles of the single particle Green's function G(uj, k) | |58| |. In the absence of any 
phase transition, G(u, k) has a form very similar to non-interacting Green's function. One can 



write [58] 



G(u,k) 



— + $(uj,k) 



(2.23) 



lo - e k + iT k 

where Z k and are quasiparticle residue and quasiparticle lifetime respectively. The first term 
in the right hand side of ( |2.23 ) is similar to the non- interacting electron Green's function 

1 



G^(uj,k) 



(2.24) 
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where 5^ = 5 sgn(e& — (j,) with < 5 -C 1. The second term, &(u>, k), represents the incoherent 
part of the Green's function. For most calculations, the first term in ( 2.23| ) is sufficient to give 
correct physical behavior up to the desired accuracy. This explains why Fermi liquid theory, 
which is very similar to non-interacting electron theory, works so well for complex systems like 
metals. In any calculation however, one has to make sure that the incoherent part of the Green's 
function does not make significant contribution. 

In terms of quasiparticle self-energy T,^(u>), and tj, are give by 68] 



jT 1 = -Z k ImS^(e fc ) (2.25) 
d_ 
doj 



z k=[l- — ReS^r(w) ) (2.26) 



We can see from ( 2.24 ) that the quasiparticle lifetime in non- interacting systems is infinite 
(rfc — ► oo). In interacting systems however, is finite. In fact, one can show that [58] for a 
3-d Fermi liquid r^T 1 oc (e^ — /x) 2 , and therefore the lifetime diverges as the particle gets close 
the Fermi surface. The closer to the Fermi surface, the more stable the quasiparticle state 
is. For states with momenta far away from the Fermi surface, the quasiparticle concept is not 
meaningful as emphasized before. 

For small u, the self energy behaves like 

ImSp^ (uj) ~ uj 2 and Re£ Pj? (cj) ~ ui (2.27) 

Substituting into ( 2.26| ) we find Z (= Z^ =p ^) a nonzero positive constant (always < 1). Z 
measures the amount of the spectral weight accumulated in the quasiparticle peak pM] . It has 



other meanings as well. It can be shown |6£, |5^] that the occupation number is always 
discontinuous at the Fermi surface with a jump exactly equal to Z. Thus the Fermi surface 
concept survives even after turning on the interaction. The existence of the Fermi surface, in 
fact, is equivalent to the validity of Fermi liquid theory. When Z = 0, as it is in 1-d Luttinger 
liquids, the occupation number doesn't have a jump at p F . Consequently, there is no well 
defined Fermi surface; although higher order derivatives of still are singular right at the 
position of the Fermi surface. Fermi liquid theory therefore breaks down in such cases. 
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The dressed interaction between two quasiparticles is given by 4-point vertex function 
r( 2 )(Pi, P 2 ; P 3 , P4) which is related to the 2-particle Green's function by 

G^(P l ,P 2 ;P 3 ,P i ) = (2tt) s G(P 1 )G(P 2 ) [«5(Pi - P 3 )S(P 2 - Pa) - <J(Pi - P 4 )<5(P 2 - P3)] 
+ i{2n) 4 5{P 1 + P 2 -P 3 - Pa) G(P 1 )G(P 2 )G(P 3 )G(P a ) T^(P 1 ,P 2 ;P 3 ,Pa) (2.28) 

with P = (v,p). Since in Fermi liquid theory, interaction between two quasiparticles is described 
by fppi, it is natural to expect connection between f pp r and T^ 2 \ To see the connection, let us 
define V(P,P'-K) by 

r(P, P';K) = (P, P';P + K,P' - K) (2.29) 

where K = (uj,k). It has been shown that the following integral equation holds under certain 
assumptions |30| 

r(p,p' ; iO = r^(p,p') + |^ / dnr^(p,Q) q ' k ^ t(q,p';K) (2.30) 

{^r J uj- q-k 

where Z is the quasiparticle residue, Q = (fJ,,p F q), with \i being the chemical potential, and 
the function is defined by 

r M ^ p/j = Hm lim p/^ ^.31) 

The order of the limits plays important role here (we will see this shortly). The integration in 
( 2.30| ) is over all directions of q. The function is closely related to the Landau's /-function 
|3(|. More precisely, at low energies putting all the momenta on the Fermi surface we have 

f pp , = Z 2 T^(p,p') (2.32) 

The notation means that P = (fJ-,p F p) and P' = (fj,,p F p'). f pp i does not have direct physical 
meaning by itself fl6q| . However, one can change the order of limits in ( 2.31| ) to define 



r (fc) ( P) pi) = lim lim r (p p' j (2.33) 

Z 2 T^ is actually the physical forward scattering interaction between two quasiparticles p6| . 
Taking the limit k — > and — > in ( 2.30| ), one can find the relation between Z 2 T^ and 
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fpp' 

„2 



Z 2 F^(p,p') = f pp , -J^JdQ fpq Z 2 T^{q,p>) (2.34) 

The quantity Z 2 T^ is usually called scattering amplitude or T-matrix [56, 43]. More precisely 
the T-matrix t pp '{q,uj) is defined by 

t pp ,(q, to) = Z p Z p ,T(P, P'; Q) (2.35) 

A modified version of ( 2.30|) gives t pp /(q,u>) in terms of f pp > [56] 



t pp/ (q,uj) = + V fpk .} Vk , ri t kp >(q,uj). (2.36) 
~ q ■ v k — u> + to k 

k 

This equation is usually known as Bethe-Salpeter equation. One can again write 

t a € = + 4a • a'ttj (2.37) 

pp' pp' pp' v ' 



In the limit q —* 0, u/q —* 0, and for p and p' on the Fermi surface, Eq. (2. 36]) gives 



f- Fl 



l + F l /{2l + l) 

tf = i + z x i{% + 1) (2 - 38) 

where if and tf are spherical harmonics. 
2.3 Fermi Liquid Theory in 2-d 

In the last section, we introduced Landau's Fermi liquid theory as a powerful tool to describe 
fermionic systems in 3-d. It is well known, on the other hand, that 1-d electron systems are not 
Fermi liquids but Luttinger liquids (7(J. The question for 2-d systems is more subtle and has 
been a subject of controversy for almost a decade. Anderson in 1990 [ p7| , |38| | suggested that 2-d 
interacting electron systems are also Luttinger liquids and not Fermi liquids. He conjectured 
that Landau's quasiparticle interaction function would become singular in 2-d. The singularity 
he suggested had the form 

p- (p-p') 

Up' — rz — ( 2 - 39 > 

\P~ P\ 
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Figure 2.3: Crossed channel diagrams contributing in I pp i. 

This singular interaction will result in vanishing quasiparticle residue and break down of Fermi 
liquid theory |43|| . One immediately notices that the numerator of ( 2.39| ) vanishes (as p — > p') 



when both particles are located on the Fermi surface. Thus the singularity can be detected only 
when the particles are away from the Fermi surface. Methods such as Renormalization Group, 
which put their initial momenta on the Fermi surface, are unable to detect this singularity. 
Since Anderson, there have been several attempts to prove or disprove Fermi liquid theory in 2- 
d H, 0, 0, H, [H ; all of them in agreement with Fermi liquid theory (at least for systems with 
isotropic Fermi surface). On the other hand, it is well known that Van Hove singularities |73|] 
or nested Fermi surfaces |73, 44] do not preserve Fermi liquid theory. Therefore the geometry 



of the Fermi surface plays an important role in the Fermi liquid behavior of the system. In 
this section, we shall study the effect of local geometry of the Fermi surface on the interaction 
function f pp i. We especially find a singularity at points where the Fermi surface has zero 
curvature. This singularity is shown, in the next section, to change the Fermi liquid behavior 
of the quasiparticles near the inflection points. 



2.3.1 Calculation of f pp / for a Curved Fermi Surface 



Let us start with a 2-d Hubbard Hamiltonian 



H = E e P°U^ + ^ E^ 4+fr c I?- g - i c p'! c PT ( 2 - 4 °) 

IT T* -T* ..I 
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where defines the "bare" band structure and O is the area of the sample. We wish to find 
fpp'' U P ^° secon d order in small coupling U, or more precisely, the renormalized coupling U, 
the same way as we did in Sec. (2.2.1). Note that by small coupling we mean U/t< 1 (this 
is not necessarily a low-density limit, but it is a weak-interaction limit). Following the method 
we described in Sec. (2.2.1), we find 



VP 



(h 2 g pp/ /m)5 a - a > - (h 2 gpp</m) 2 (I pp/ + C pp , ) 



(2.41) 



where the dimensionless coupling g pp > is given by [45, 43 1 

mU m 



9pp' 



ti 2 



h 1 



T7 2 



1 



p'+q P-i 



(2.42) 



we also identify the "Cooper" channel and "crossed" channel by (cf. Fig. 2.3 ) 



Ipp> 



E 

k 



p' k 



(-P ~ £J + e t 



p+p' — fe / 

h (p <-> p') 



(2.43) 
(2.44) 



-p' ' k k+Q / 

The Cooper term, C™, , is singular when p = (p + p')/2 — ► 0, but regular near Q = p — pf = 
(as is g pp '). We thus drop Cw' and approximate g pp ' « 50 • ^pp'> on the other hand, may become 
singular when Q — > 0. This can be seen by the fact that the denominators in ( |2.44j ) vanish at 
Q = 0. We therefore have to examine I pp i more carefully. Let us define 



v t =^% and (M^ 1 ) i j = 

k dk K k Jj 

Expanding the denominators to second order in Q, we have 
d 2 k 



dkidkj 



(2.45) 



Ipp' 



Til 



llT 



(2nf 



1 



v k )-Q-\Q- M- t l ■ Q (v p -v k )-Q + \Q- Mz l ■ Q 



(2vr) 2 



d 2 k- - oc UQ) 



[(v p - v k ) • Q} 2 - A 



k,Q 



where 



HQ) 



d 2 k 



[(vp-v k )-Q} 2 -Al Q 



(2.46) 



(2.47) 
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Figure 2.4: The two situations considered in the text. In (a) we assume the Fermi surface has finite 
curvature at k = k ; in (b) the curvature is zero when k = fc (inflection point). Hatched lines show the 
integration regions discussed in the text. 



The " gap" A^ ^ is denned by 



Q-M^Q 

oc \Q\ 



(2.48) 



k ' Q 2\Q\ 

and Q = Q/\Q\. For a quadratic dispersion relation = k 2 /2m, -M^ 1 (= r/m) is independent 
of k. In more general cases however, it has fe-dependence but this does not affect our calculation 
(we will consider the fe-dependence of M^ 1 when appropriate later). Thus for the moment, we 
can safely neglect the ^-dependence and write Ag^ = A^. 

Our goal is now to understand the behavior of Ig(Q) as Q — > 0. Infrared singularity occurs 
when (vp — Vk) • Q = 0. This condition is usually satisfied on a line in /c-space inside the Fermi 
surface (corresponding to a straight line in u-space). Let's confine our integration to a region 
within a cutoff A around this line, where \Q\ <C A <C p F , and p F is the Fermi momentum. The 
integration inside the Fermi surface then has the form 

' A du 1 JA-A;jl i 2A. 



A U z 



Q 



1 |A - A c 

a-, H WTa c 



Q • 



A, 



A 



2 
A 



(2.49) 



Thus the integration inside the Fermi surface contains no singularity as A^ — ► 0. Singularity, 
if any, should therefore come from the integration near the Fermi surface. Let ko be a point on 
the Fermi surface, satisfying (vp — Vk ) ■ Q = 0. We first consider points where Fermi surface 
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has finite curvature (Fig. |2.4| (a)). Writing k = ko + k±e± + k\\i.\\, where en (e±) are unit vectors 
parallel (perpendicular) to the Fermi surface at ko, we then define eg by the expansion 



e£ = E F + Vk kj_ + ak\ + bkj + ck±k\\ 



Near ko therefore, can be written as 



Vk = v ko e± + 7||fc|| + j±kj_ 



where 



(2.50) 



(2.51) 



7_i_ = 2ae± + cen 



7|| = c e± + 26eii . 



(2.52) 



Fermi surface also is defined by 



ffc ^± + + + ck±k\\ = 



(2.53) 



Since the Fermi surface has finite curvature at ko, 7y and b are non-zero. Rescaling the integral 
in ( 2.46| ) so that Vk = b = 1, we have a Fermi surface near ko obeying ~ — k± (cf. ( |2.53| )). 
It is not difficult to see that singularities in I pp > , for p, p 1 close to ko , arise only when 



7|| • Q - 



In that case: 



A 2 



dk\\ 



7T 



\Q\^ 2 



(2.54) 



(2.55) 



To get ( 2.55 ) we have made use of the following integral 



du 



u 



— i / \ 

n 2 -A2 VA^^^ + VAI +tan VA j 



:(-ln- 



(2.56) 



This singularity is just that found already for a circular Fermi surface [42, 43]. It only occurs 
if p and p' approach each other along the direction orthogonal to 711. For a circular Fermi 
surface, this direction is orthogonal to the Fermi surface. For non-circular Fermi surfaces (with 
no inflection point) on the other hand, 711 -Q = for directions which are not necessarily normal 



Chapter 2. Fermi Liquid Theory 



23 



to the Fermi surface; but the singularity is still along lines. This case was also studied away 



from the Fermi surface by Fukuyama and Ogata [72]. A closer look at the singular region in 
fc-space for the case of a circular Fermi surface shows that it is bounded by circular arcs of 
opposite curvature which touch at the Fermi surface E3]. 



This singularity does not change the Fermi liquid behavior 43]. A simple way to see this, 
is to look at the quasiparticle energy 

tp = 4+Y,fpp> dn p ( 2 - 57 ) 

Singularity in f pp i is summed over in this equation. In our case here, the [Q] 1 / 2 singularity is 
always along a line. Integration over this term removes the singularity and gives Fermi liquid 



behavior as one approaches the Fermi surface [42, 43]. 



2.3.2 Calculation of f pp / near an Inflection Point 

At points where the Fermi surface has zero curvature, 711 — ► and consequently, condition 



(2.54) is always satisfied. It is therefore natural to expect a singularity in any direction of Q; 
which may now cause break down of Fermi liquid theory. We wish to study this singularity 
near an inflection point. We must therefore extend (|2.53| ) to higher order, and also include 
the /c-dependence of (which now vanishes in certain directions). Expanding about an 

inflection point at k = ko, we write 



g = v ko k ± + aki + gkS (2.58) 



(cf. Fig. 1(b)). Substituting directly into Q2.44 ) we get 





i—gkpv ko 1 


/>'. 




'-gAyv k0 \ 



h(Q) = 77^ / d h / ' " I rr^; t^z I + (Q - ~Q) ( 2 - 59 ) 



We have defined k = p—ko (i.e., the center of mass momentum relative to ko), c Q = 2aQ±/3gQ\\ 
and 

Ql 0, 

T(k, Q) = k\ - -JL + c Q (« ± + Z±) (2.60) 
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Integration over k± is straightforward. One finds 



ip(Q) 



l 



2aQj 



dku In 



—A 



fe 2 + Qu fey - (c Q g/v ko )kf { + T(k, Q) 



+ (Q - -Q) 



(2.61) 



kf + Qfy - {c Q g/v ko )A 3 + T(k,Q) 

The fey term in the numerator is small compared to the other terms unless if Qy — > (or 
Cq — > oo). We therefore neglect this term except when Qu <C Q± (we will come back to this 
case later). To calculate the above integrals we should use 



J dum\u 2 + (3u + 7I = (u + /?/2)ln|n 2 + /3u + 7] - 2u 

e { -w)\ m f" + ^ 2 + ^ + W tan- («+££ 



+2J\W\ 



(2.62) 



where = 7 — /3 2 /4. However , it turns out that all the terms in the above integral are regular 
as Q — > except for the "tan -1 " term. Keeping only the singular terms we find 



W 



+ (Q 



where 



Qi 



Wtf = T(H, Q) + = — k 2 — c Q ( K± + Q ± /2) 
In the limit where k, Q <C A we approximately get 

7T 



- [| W + | 1/2 ^(W_4_) - |^| 1/2 6>(W-)] + ... 



-Q) 

(2.63) 
(2.64) 

(2.65) 



where (...) represents non-singular terms and W± = W^^a. 

This behavior of I pp i is rather complicated. We can get some feeling for it if we first take 
the case where p = ko, i.e., where one quasiparticle is at the inflection point (see Fig. |2.5[ ), 
Then for most regions in fe-space, ( |2.65 ) gives 

hoM+Q \Q^T2 ( 2 - 66 ) 

This is the same inverse square root as found for the circular Fermi surface (cf. ( 2.55] )), but 



now extending over large "butterfly wing" regions in fe-space.) Fig. 2.5 shows different regions 
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Figure 2.5: Three different regions in fc-space around the inflection points with different singular behav- 
ior. We have put p = k therefore the regions shown are actually the regions where p' can be. Between 
I and II, there is a crossover whereas between II and III there is a sharp edge. 

in fc-space with different singular behavior. The singularity in the left wing of the butterfly 
region comes from switching pand p', i.e. putting p' at the inflection point. 

When Qii — ► 0, the above approximation does not hold anymore. In that case ( |2.5S ) becomes 



ip(Q) 



l 



2aQ± J-a 

CK/fl) 1/3 

2aQ± 



dk i 



K- A/9) 1/3 

-K>fc±/s) 1/3 



1 



k± - p± + k ± 



(k ± + Q ± /2) 1 / 3 -(k ± -Q ± /2) 1 



/3 



+ (p <-> p' ) 



+ 



(2.67) 



where ( dzj (z 3 — 1) = — it/ y/3. If we again put p = ko, we get 



ip(Q) ~ Q[ 



-2/3 



(Qll < Qj 



-1/2, 



(2.68) 

There is a crossover 



The singularity here is stronger than what we found before (~ Q 
between these two different behaviors (cf. Fig. |2.5| ). 

It is very important to note also the non-singular region (region III in Fig. |2.5| ) which 
lies close to the Fermi surface as Q approaches the origin. There is no singularity in the 
interaction between 2 quasiparticles if they are both on the Fermi surface. Thus any calculation 
(perturbative or renormalization group) of the scattering between quasiparticles confined to the 
Fermi surface will miss all the singular behavior in ( p. 66 ). 

In the more general case where neither particle is at ko, the behavior of /„„/ is not so easily 
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displayed. As the quasiparticles approach each other one finds the simple behavior 
2.4 Calculation of Self Energy 

Fermi or non-Fermi liquid nature of the system can be checked by calculating the self energy 
We already have seen how the quasiparticle lifetime and quasiparticle residue are related 
to the real and imaginery parts of the self energy (cf. ( |2.25| ) and ( 2.26[ )). The imaginary part 



of T,fi(uj) is specially important because it can be calculated easily using the Fermi golden rule. 
We saw that in a 3-d Fermi liquid, at low energies, we have 

ImSp F (u) ~ u? (2.70) 

In this section we first see how this relation is different for a 2-d Fermi liquid with non-singular 
interaction. We then find the effect of our singularities on ImS Pf (w) and show the deviation 
from Fermi liquid behavior. 

2.4.1 Non-Singular Interaction 

Let us first consider the non-singular case. Landau's /-function is therefore regular and given 
by the cylindrical harmonics // the same way as before. For the present calculation, it is enough 
to consider the isotropic harmonic and write 

f PP > = fo (2.71) 

Having Landau's /-function, we can calculate the quasiparticle T-matrix using the Bethe- 
Salpeter equation ( |2,36| ). In the isotropic channel, which we consider here, the T-matrix is 
simply given by [Eg] 



t pp/ (q, to) = t ( V ) = (2.72) 



A 

l + /ox°(<Z> 

where yp(q,uj) is the 2-dimensional Lindhard function 



X °(q,u;) = N(0) 



1 nf ^ 2 - 1 ) | A 1 -^ 



(2.73) 
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r] = Lo/qv F and density of states at the Fermi surface iV(0) = p F /2nv F . All the physical quan- 
tities can then be found from the T-matrix. What we wish to do here however, is to calculate 
the quasiparticle self energy and compare it with the 3-d case discussed before (calculation we 
present here is based on Ref. |43[| ). 

It is always easier to calculate the imaginary part of the self energy which is equivalent to 
the scattering rate and is given by Fermi golden rule [56] 

We calculate the energy on shell (i.e. u = ep). Assuming q <C p F we have 

tp-q = ep- v F qp = u) - v F qp (2.75) 
where p = p ■ q. The 5-function in ( |2.74| ) can be written as 

/oo 
di/£(i/ - u + ep-$)5(u - ep - e p -, + -)) 



CO 

oo 



dv5{y — qv F p)5{v — qv F p) (2.76) 
with fjf = p' ■ q = cos 9'. On the other hand (at T = 0) 

nj(l - n°, +? .) = 0(p F -p')%' + — Pf) (2-77) 



Eq. ( 2.77 ) requires p' < p F and fjf > 0. Therefore 

2^ n ,? 1 - n g+a> = ~o — = tt~\2 / ^ de ( 2 - 78 ) 

p J- n/2 2tt J p qiM > 2tt (2tt) 2 J- n/2 



Substituting all back into (2.74) we find 



Im£» = W I = l f (g» " <P> F IJ)0{?Or ~ ~ QV F ix) (2.79) 

4TTV„ J U.2 a 2 _ j, 2 



One can write 
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Because of the 5- function, the //-integral will require /i = v /qv F = ij. As a result 

bnEpCw) = [ qdq [ ud J^!'^ 0(lo - u)9(qv F - u) (2.81) 

\2-k) 6 v f J J q z v p — v A 

In the isotropic case, t(q, v) is only a function of r\ (cf. ( |2,72| )). For rj < 1, which is actually 
the case here, we have 

f / \ fo _^ A) 

W 1 + / (1 - VvT^)JV(O) iV(0) 1 + iAvily/l^ 1 ' j 



where A = -Fo/(l + -Pb) with i<o = N(0)fo dimensionless Landau parameter. In this case ( 2.81 ) 
becomes 

>' f^dut (2.83) 



(2irv F ) 3 Jo J v / Uc rj 1 - rf 
where oj c = Av F is the energy cutoff. Substituting t(rj) from ( |2.82j ) we get 

Im£ p (w) = K vdv f 1 - —r^ 7s 2 ^K T udu f (2.84) 

where iT = j4g/27ru F p F . In the last step we have dropped the tj 2 term in the denominator. We 
therefore simply get 

ImEp(w) ~ / dvvln\ — \ ~ w 2 (l + 2 In I — I) (2.85) 
Jo v UJ 

The leading u dependence is therefore 

ImS p (o;) ~ w 2 lnw (2.86) 



Although this is different from the 3-d case (2.70), it does not change Fermi liquid behavior 



[43]. Now let us see how our singular interaction will change this result. 



2.4.2 Singular Interaction 

We now calculate the self energy for the singular interaction we found near inflection points. 



As in some previous studies of various kinds of singular interaction [03, 42, 4C, 75, 76], it is 
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Figure 2.6: Self energy to second order in A p 



p-q- 



convenient to begin with that part of T,Juj) which is 2nd-order in the singular interaction. Let 



us write 

Up' = fp P ^ T + \jp' (2.87) 

where f pp ^ T is the non-singular part of f pp i and X pp i is the singular part. We calculate S^w) to 
second order in X pp i, but to all orders in f p J" T - Since we are not particularly interested in the 
angular dependence of the regular Fermi liquid contribution, we will simply let f pp ^ T -> /o, a 
constant as before. We do our calculation for p = Icq, i.e. at the inflection point. The singular 
term X pp i therefore is given by (|2~66; ) and Q ). Since ( ^68|) is only valid in a region very close 



to the normal axis, we take X pp i every were to be 

1 



Xpp> ~ , (2.88) 



We again use ( 2.74| ) as our starting point (note that ( p.74 ) is valid for u > 0, as can be 



immidiately seen by the form of the occupation numbers. A full version is given in Ref. |36j 
In order to calculate to 2nd order in X pp ', \t pp /\ 2 in ( |2.74 ) should be replaced by [^] 



A 

\tpp>{q^)\ 2 ~ ^ (2.89) 

where x°(Q;^) is the generalization of ( |2.73| ) for a non-circular Fermi surface. This is equivalent 
to calculating the diagram presented in Fig. ^1] [43], where xi.Qj 1 ') 1S the usual fluctuation 
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propagator 

Since we are going to calculate everything near the inflection point, we can assume that x°(<f, co) 
behaves regularly at the inflection point and treat it as a constant. At T = we have 

We expand near the inflection point ko by 

efc = Vfco«± + »«| (2-92) 
where k = k — k^. From now on, for simplisity, we take Vk = g = l.Thus 

ep = -Q± - Qf\ 

£p-q = ~<1± - Q\\ 

zp + f=q±-Q± + (q\\-Q\\) 3 (2-93) 

The self energy now becomes 

Im£,,» ~ / / d 1\\ / A Q3 d Q± d 1±^ -Q±-Qf\+Q±+ Q\\) 

x6(u-3Q$q { \+3Q\\4) (2.94) 
Performing the q± integral, we get 

ImS fc0 (a;) ~ J ^ J dq\\ £ ^ dQ ± (u-Q ± - Q|)0(w -Q±- Qjj) 



X( 



5{u-3Qfq\\+3Q\\<4) ( 2 - 95 ) 



Integrating over Q± is also trivial 

v dQ\\ r\ 



Im£ fc » ~ j ^ j^dq 6(u> :\Q-<, ■ SQ qf\) (2-96) 

We now change the variables to Q\\ — > w 1 / 3 a; and gy — > uj l ^y. Integration limits also change to 
ztA/w 1 / 3 . For small uj (<S A 3 , or more precisely <C gA 3 ), the integration limits can be replaced 
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by ±00. However, one has to make sure that the integral does not diverge as the limits go to 
infinity. We finally get 

poo J T poo 

ImE fco (w) ~ w 4 / 3 / — / dyd(l - Zx 2 y + 3xy 2 ) (2.97) 

J— 00 J— 00 

The integral is regular at integration limits and also at x = 0. Therefore we can confidently 
write 

ImS fc[) (w)~w 4/3 (2.98) 

This is clearly different from the Fermi liquid form ( |2.70|) or the 2-d version of it ( |2.86|) . This is 
also different from the result obtianed by Fukuyama and Ogata W% which they find the power 
to be (5/3) instead of (4/3) for inflection points. The reason for this descripency is that they 
miss the singular interaction in their calculation. Notice that if we drop the 1/Qn factor (which 
comes from the singular interaction) in the integrand of ( |2.96| ) , we get ImS ~ cj 5 / 3 in agreement 



with Ref. [72|. 



We can repeat this calculation for a curved Fermi surface using 

e k = K± + K jj (2.99) 



instead of ( |2~9p . Eq. ( p^ ) then becomes 

cA r A 



/A rj\ 
d Q\\ dq»S(u + 2Quq» - 2qh 
-A J-A 



-A J-A 

A dQ\\ 



-A \Q\\\ 

uj 2 Itluj (2.100) 



which agrees with (|]8|). Although (^9^ ) is different from ( p.l00|) , the difference is not enough 



to break down Fermi liquid theory (in the sense that Z / and Fermi surface exists) |f72j| . On 
the other hand, higher order calculations can change this result significantly. For example simple 
estimate indicates that the 4th order calculation gives marginal Fermi liquid behavior ImS ~ u. 
Systematic higher order calculation therefore is necessary to acheive the final conclusion Jnj] . 
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2.5 Summary 

Our calculation of f pp i for a 2-d electron system reveals a 1-d singularity near curved points 
of the Fermi surface as p — > p' . This is in agreement with previous results and is shown to 
preserve Fermi liquid behavior Near inflection points of the Fermi surface, on the other 
hand, we find a 2-d singularity again with \p — p / \~ 1 / 2 behavior. To second order in the singular 
term, the imaginary part of the self energy shows an w 4 / 3 -dependence in disagreement with the 
ordinary Fermi liquid result (i.e. u; 2 lnu;). This however does not invalidate Fermi liquid theory 
f72f . Inclusion of Higher order diagrams is necessary for any rigorous conclusion. 



We should point out that our result is different from Ref. \72\ in a sense that they do not 
incorporate singular interaction at the inflection points in their calculation of the self energy. 
They report u 5 / 3 -dependence for the imaginary self energy near an inflection point, different 
from what we find (u 4 / 3 ). Ref. [78| has also investigated the effect of the inflection points. They 
however consider the scattering between two quasiparticles located on two inflection points 
opposite to each other with respect to the center of the Brillouin zone. This is clearly different 
from our calculation. Employing renormalization group technique they find that inflection 
points would increase the lower critical dimension for instability from Fermi liquid fixed point 
from d=l to d=3/2. Thus 2-d systems will still be stable under such an effect. 



Chapter 3 



Nearly Antiferromagnetic Fermi Liquid Theory 



At zero doping, all cuprates are Mott insulators with antiferromagnetically ordered spins. On 
the other hand, the physical behavior of overdoped samples is well described by Fermi liquid 
theory. However, the most interesting region in the phase space of high T c 's is none of the above 
phases, but the intermediate doping regime. One can think about the intermediate doping as a 
continuation of the metallic phase or the extension of the insulating phase. Depending on which 
way to choose, one may end up with two different schools of thought based on Fermi liquid 
or non-Fermi liquid theories. At this time not only no consensus about the older theories of 
high T c has been achieved, but also different new theories such as, SO(5) theory jnj, theory of 
stripes [80], and more recently theory of nodal liquids [^lj, etc. have emerged. In this chapter 
we shall focus on one of the older theories which is based on the Fermi liquid assumption, and 
because of its simplicity in calculations, has attracted considerable attention especially among 
experimentalists. 

Nearly antiferromagnetic Fermi liquid theory (NAFL) is a phenomenological theory moti- 
vated by nearly ferromagnetic spin fluctuation theory or paramagnon theory, proposed in 1960's 
to describe metals near their ferromagnetic instability as well as liquid 3 He. Both theories are 
based on the assumption that the excitations of the electronic system are Fermi liquid quasipar- 
ticles, interacting among each other via the damped spin waves. Besides the similarities, there 
are fundamental differences between these two theories which have made NAFL not as widely 
accepted as paramagnon theory. In this chapter, after a brief review of paramagnon theory, we 
introduce NAFL emphasising on the similarities and differences between the two theories. 
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3.1 Paramagnon Theory 

Fermi liquid theory is primarily constructed to describe properties of low energy electron-hole 
excitations |p8[| . There also exist other excitations in a Fermi liquid, which are inherently 
different from electron-hole excitations (although they can be views as collective resonances of 



them [68]), such as density oscillations and damped spin waves. These excitations are usually 
known as plasmons and paramagnons respectively. On the other hand, although thermodynamic 
and static measurements such as specific heat and spin susceptibility can be expressed in terms 
of the Fermi liquid parameters, the information contained in Landau parameters is not enough 
to describe some other measurements, e.g. transport properties and superconductivity (or 
superfluidity in 3 He). In particular, the interaction between the quasiparticles plays important 
role for these properties. We saw in the last chapter that the effective interaction between 
quasiparticles can be obtained from the Landau interaction function via Bethe-Salpeter equation 
( |2.36j ). This is however, not always the easiest route to take. The other way is to assume a 
model Hamiltonian for interaction between quasiparticles and apply field theory techniques to 
calculate the physical properties. The reliability of these calculations should then be judged by 
experiments. One of the theories that follow this path is paramagnon theory which we describe 
briefly here. 

Magnetic susceptibility of a Fermi liquid is given by |5Q] 

X - = 4vr 2 (l + Zo) (3 - 101) 



If the Landau parameter Zq — * —1, the denominator of ( 3.101| ) vanishes and therefore magnetic 



susceptibility diverges. The system in that case undergoes a phase transition to a ferromagnet- 
ically ordered state. When the denominator in ( [3.101 ) is not zero but very small, the magnetic 



susceptibility is still large and the neighboring spins still have the tendency to align, although 
the long range order no longer exists. On the other hand, the interaction between the quasi- 
particles is given by the T- matrix which can be found by Bethe-Salpeter equation ( p. 36 ). In 
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particular ( [2.38 ) gives 



tg = V(i + z ) 



(3.102) 



which also diverges as Zq — > — 1. Thus strong magnetic susceptibility results in strong in- 
teraction between quasiparticles; which tends to align the neighboring spins. The dynamic 



(3.103) 



susceptibility of a Fermi liquid was suggested by Leggett [55] to be 

X[q ' UJ) l + Z oX (°)(q,u)/N(0) 
where 

x (o) is 

the susceptibility of the non-interacting electron system (particle- hole bubble). In 
3-d isotropic systems, 1S given by Lindhard function [68]. 

Eq. ( |3. 103 ) was proposed by Leggett using deductive reasoning. An attempt to derive it 



was first made by paramagnon theory [82]. The theory assumes that the (effective) interaction 
occurs only between quasiparticles with the opposite spins. The model Hamiltonian for this 



interaction was first proposed by Berk and Schrieffer [83] as 



Hint = I d r n^(r)ni(r) 



I 

V 



Eel ^crct cti 
k+q * k'—q & 



(3.104) 



k,k',q 



where cL (eg) is an operator that creates (annihilates) a quasiparticle with momentum k and 
I > is the coupling constant. Performing Eliashberg type calculations, Berk and Schrieffer 
showed that in nearly ferromagnetic metals, such as Pd, superconductivity is suppressed as a 
result of ferromagnetically enhenced coupling between the quasiparticles. This model was then 
adopted by Doniach and Engelsberg f82[| and Rice for liquid 3 He. 

Spin susceptibility can be calculated from ( 3.104] ) using Kubo formula [68] and RPA ap- 
proximation 



Comparing (13.1030 and ( |3.105| ) we find 



x (0) (g» 

1-I X (°)(?» 



(3.105) 



Z = -N(0)I 



(3.106) 
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Thus the interaction I is closely related to Fermi liquid interaction function f pp /. In low energy 
and momentum, one can expand ( |3.105| ) to get 

w) = i-JTWW^WJW^ (3 - 107) 

where I = IN(Q) and q = q/p F . The susceptibility x(^<^) m (|3.107 ) is peaked at q = 0, 
reflecting the proximity to ferromagnetism. Near the antiferromagnetic instability, one expects 
the susceptibility to be peaked near q = (±7r/o, ±7r/a, ±7r/a) where a is the lattice spacing. 

x(q, to) can be viewed as the propagator of a bosonic field that mediates the interaction 
between the fermionic quasiparticles. The corresponding bosons are called paramagnons which 
are critically damped spin waves. Eq. ( |3.107j ) can be written as 

X (q, u) = f , (3.108) 

1 + 4T<T - tu/oJ q 



where xo = SN(0) is the enhanced static susceptibility, £ = ySI/12p 2 F is the ferromagnetic 
correlation length and co q = Aq/nIp F is an energy scale that characterizes the damping of the 
bosonic particles. Here, S = 1/(1 — /) is the Stoner enhancement factor which measures the 
proximity to the ferromagnetic instability. 

Although paramagnon theory gives the right form for the spin susceptibility, it has some 



internal inconsistencies [ S3 1 . For example if we calculate the density-density correlation, instead 
of ( plPCD we get g|] 

F = N(0)I (3.109) 

( 3.106 ) and ( |3.109 ) can be satisfied simultaneously only if Fq = —Zq. This is obviously not 
always the case. For example for liquid 3 He [58], Fq = 10 but Zq = —0.67. Thus paramagnon 
theory cannot be the right quantitative theory for the bare quasiparticles in liquid 3 He. Even 
without comparison with experiments the theory has internal inconsistencies. For example RPA 
approximation does not satisfy Pauli principle |54j . Another problem is the absence of Migdal's 
theorem in the electron-spin fluctuation interaction, first raised by Hertz et al. |85|| . Vertex 
corrections was shown to be of the same order as the bare interaction. This invalidates the 
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Migdals theorem and therefore any Eliashberg type calculation. We will study this point in the 
context of NAFL theory later. Besides quantitative inconsistency, paramagnon theory gives 
some qualitatively right results. We will not discuss the applications of paramagnon theory in 
this thesis. Instead, we proceed to the antiferromagnetic version of the theory which is relevant 



to high T c superconductivity. Interested reader should refer to Ref. [63] for more information 
about paramagnon theory. 

3.2 Antiferromagnetic Spin Fluctuations 

It is more or less well established that high T c materials are effectively 2-dimensional electron 
systems close to the antiferromagnetic instability. As we showed in the last chapter, a two- 
dimensional electron system with a smoothly curved Fermi surface can be a Fermi liquid (at 
least in weak interaction regime). On the other hand, NMR experiments [jl4| and inelastic 
neutron scattering (INS) |8(| indicate enhanced staggered susceptibility, even at optimal doping. 
One therefore might be tempted to generalize paramagnon theory to antiferromagnetic systems 
to describe high T c superconductors. It is well known that for a nested Fermi surface, the RPA 
susceptibility (3.105) diverges at the wave vector q = Q = (±7r/2a, ±7r/2a) resulting in spin 



density wave instability [44]. For a general Fermi surface however, microscopic derivation of 
the spin susceptibility fails to succeed because the incoherent parts of the quasiparticle spectral 
function play important role in the calculation of the real part of the particle- hole bubble [87]. 



The low frequency part of the imaginary part of the particle-hole bubble, on the other hand, is 
given merely by the coherent parts and can be easily calculated. We will discuss this issue when 
we talk about vertex corrections later. Based on NMR and NIS results, a phenomenological 
susceptibility analogous to (|3.10§| ) was suggested by Millis, Monien and Pines [23] for high T c 
systems 

X(q, w) = _ *Q , - (3.110) 
l + £ 2 (q-Q) 2 -iu/u s f 




Figure 3.7: Fermi surface in the first Brillouin zone, with the value for t and if given in the text; and 
we assume n=0.75. Calculations are presented here for the wave vectors k\ and kh- 

where xq ^ Xo is now the static staggered susceptibility, £ is the antiferromagnetic correlation 
length and u s f specifies low frequency relaxation behavior. As in paramagnon theory, quasipar- 
ticles are again assumed to interact with each other via the spin fluctuations. The Hamiltonian 



that describes this interaction was written as [31 1 



2 k+q,a 
P.°" q,k,a,/3 



where g is spin-fermion coupling constant and 

ep = — 2t[cos(p x a) + cos(p y a)] — At cos(p x a) cos(p y a) — \i (3.112) 



is the quasiparticle dispersion relation [35, 88|, with t = 0.25eV and if = — 0.45t. S(r) is the 



spin fluctuation density operator with a bosonic propagator given by ( p.llO| ). 

The presence of if allows "hot spots" pj| on the Fermi surface (Fig. |3.7| ) which are points 
that can be connected by Q; this leads to singular behavior when £ in ( 3.110Q diverges p9| , 



j9|. Quasiparticles near the hot spots are highly scattered and play important role in NAFL 
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calculations. Full analysis of the scattering rate 8£, [X]] indicates linear T dependence at the 
hot spots, while T or T 2 dependence at other points, depending on the temperature and the 
distance from the hot spots. This is used to explain the linear T dependence of the normal state 
resistivity of cuprates and the crossover between linear T behavior to T 2 behavior 8£, 90]. Some 



other properties of the cuprates such as superconducting T c [31, 35, 36], Hall conductivity [90] 



magnetic and scaling properties [91], |92[, evolution of Fermi surface |93|, Raman scattering 
], etc., have been calculated using NAFL theory. It is not our intention in this thesis to 
describe these calculations in detail. Instead, we would like to study the self-consistency of the 
theory. We proceed by briefly discussing how d-wave superconductivity naturally arises from 
NAFL theory. We then discuss the self-consistency of NAFL and calculations made by this 
theory. 

3.3 d-wave Superconductivity 

In NAFL theory, antiferromagnetic spin fluctuations play the same role in the pairing process 
as phonons do in conventional superconductors. One can therefore apply BCS theory, replacing 
phonons with magnons. There is however a principle difference here, that is the effective 
interaction potential 

V eS = g 2 X (q) (3.113) 

is always positive (repulsive). Here g = \J?>/Ag (the factor comes from the consideration of the 
Pauli matrices). BCS gap equation for this interaction is 

A, = -g 2 I d 2 k>A k , X (k - fcO tan t fc,/2r) (3.H4) 
Since x is positive, for an s-wave gap A^, both sides of (|3.114 ) appear with different signs and 



therefore s-wave paring is not possible. However, since x(o) is peaked at q = Q, d x 2_ y 2 pairing 



is possible because in that case A^+q = — A& and this changes the negative sign in (3.114) 
to positive. Therefore d x 2_ y 2 pairing is a natural consequence of NAFL theory. Experimental 
observation of d-wave pairing is considered as one of the triumphs of NAFL theory. 



Chapter 3. Nearly Antiferromagnetic Fermi Liquid Theory 



40 



In order to get high T c , spin fermion coupling should be strong. At the beginning it was 
thought that short quasiparticle lifetime due to the strong scattering would reduce T c dramat- 



ically |pSfl . Weak coupling [|3l| and also strong coupling (Eliashberg) |35|, |3g] calculations of 
Urbana group however, have found a T c as high as 90K, even including the lifetime effect [36]. 
Weak coupling calculation in Ref. gives a T c which is well approximated by 

T c = a C^j e' l/x (3.115) 

where V ~ 0.04ey is an energy scale, 

A = i]g 2 XoN(0) (3.116) 

is the dimensionless coupling constant, and a and r/ are dimensionless constants of 0(1). Strong 
coupling calculation of T c was also done by Monthoux and Pines (MP) p5| , 36] using Eliashberg 
formalism. For t' = 0, a relation similar to ( 3 . 11 5| ) was again obtained 



0.636 ( ^ ) e -V0.402JV(0)5 (3 n7) 



where 



IT 2 



N(0) = --J2^G R (p,0) (3.118) 

7T — 

P 



is the tunneling density of state calculated from one particle retarded Green's function at T c 



[35]. It is interesting that T c is only sensitive to g and not other parameters. Fitting to the 



measured value for T c , one can determine g. For T c = 90K, one gets jw0.8 [p^[ . 
3.4 Self- Consistency of NAFL 

Besides the similarities, there is one fundamental difference between NAFL theory and para- 
magnon theory: charge properties of ferromagnets are still metallic while antiferromagnets 



are insulators. The reason behind this difference as pointed out by Anderson [95, |9l| is that 
ferromagnetism is a natural consequence of Hund's rule: exchange interaction favors spins to 
be parallel in order to keep electrons as far as possible and make the repulsion energy mini- 
mum. Antiferromagnetism on the other hand, is a result of superexchange interaction which 
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is inherently a non-perturbative effect and is a property of Mott insulators. Going out of the 
antiferromagnetically ordered state involves a metal-insulator transition which does not have 
a counterpart in ferromagnetic case (there is a slight difference between strong coupling and 
weak coupling here. In weak coupling the Neel ordering and metal-insulator transition - which 
is not of the Mott type - happen at the same time whereas in strong coupling the system is 
already an insulator above the Neel temperature). It therefore completely makes sense to talk 
about a nearly ferromagnetic Fermi liquid but one should be very careful to talk about a nearly 
antiferromagnetic Fermi liquid. On the other hand, almost none of the properties of high T c 
superconductors resembles of a Fermi liquid, even outside the antiferromagnetically ordered 
(insulating) region of the phase diagram (except for very high dopings). Thus it is hard to 
believe that the electron system turns into a Fermi liquid state from a highly non- Fermi liquid 
insulating state by slightly adding doping. 

Another problem, which is more computational than fundamental, is consideration of the 
vertex corrections (note that in pseudogap problem the strong frequency dependence of the 
vertex correction may make a difference between having or not having a pseudogap |Mj . I that 
sense vertex correction is a fundamental problem). In order to calculate T c in strong coupling 
limit, Eliashberg formalism is commonly used. A necessary condition for the applicability 



of Eliashberg formalism is Migdal's theorem [46]; which states that the corrections to the bare 
spin-fermion vertex are negligible compared to the bare vertex. In case of electron-phonon 



interaction, the vertex corrections are of order of (ujd/E f ) ~ \JmjM <C 1, where m and M are 
the masses of electron and nucleus respectively. Spin fluctuations however are different from 
phonons in three aspects. First, phonons are external excitations while spin fluctuations are 
excitations of the electronic system itself. Secondly, unlike phonons, spin fluctuations are not 
propagating modes. Thirdly, the ratio J/E F , with J being the exchange coupling responsible 
for antiferromagnetism, is of order of one, so there is no small parameter like {u)d/E f ) in 
spin-fluctuation interaction 



Early estimates reported small vertex corrections [97[. Calculations of Ref. [48] and [49] on 
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the other hand, indicated that close to the antiferromagnetic instability point, the first vertex 
correction contains a logarithmically divergent term. However, sum of a series of divergent 



diagrams results in a power law behavior for the dressed coupling [49, 5C] 

Stat ~ g{-f (3.119) 
a 

where f3 is a small positive number (we should mention that in the original publication [5 
was reported to be negative but in later publications (5 was claimed to be positive [p0[ ). Thus 
near the antiferromagnetic instability (£ — * oo), ptot become very large. In other words, vertex 
corrections tend to increase the spin-fermion coupling. 

On the other hand, a work by Schrieffer [47] showed that right outside the antiferromagnet- 
ically ordered phase, the dressed spin-fermion vertex is given by 

^oc[(?-Q) 2 + r 2 ] 1/2 (3-120) 

Thus very close to the antiferromagnetic instability (£ — ► oo), the full vertex vanishes as q — > Q. 
The vanishing of the dressed vertex is a consequence of Adler principle: in the ordered state, 
magnons are Goldstone modes and their interaction with other degrees of freedom should vanish 
at the ordering momentum to preserve the Goldstone modes to all orders in perturbation theory 
The above statement means that the vertex corrections should be equal and opposite sign 
to the bare vertex to make the dressed vertex vanish. Consequently, Migdal's theorem can not 
exist near the antiferromagnetic instability. However, superconductivity does not happen at 
very low dopings either. The relevant question therefore is the validity of Migdal's theorem 
at higher dopings, outside the ordered state where superconducting transition occurs at high 
temperatures. 

The first attempt to calculate vertex corrections at optimal doping was made by us [98]. 
Using the same parameters as used in NAFL theory, we showed that the first correction to the 



bare vertex is actually the same order as the bare vertex [ 98 1 . The sign of the correction is also 
controversial because it determines whether vertex corrections increase or suppress the effect of 
the spin fluctuation coupling. Our finding had a sign opposite to the bare vertex. The effect of 
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the correction is therefore to reduce the spin-fermion coupling, in agreement with Schrieffer's 
result ( p. 120 ). In the next section we describe our calculation in detail. We then discuss the 
effect of inclusion of the quasiparticle renormalization in our calculation. 

To conclude this section we mention one experimental consideration. The susceptibility 



(3.110) is adopted by NAFL theory only on phenomenological basis and by fitting to experi- 
ments. Consequently, NAFL theory is incapable of explaining the incommensurate peaks ob- 
served in inelastic neutron scattering (INS) experiments on La2- x Sr x Cu04 |H| and recently on 
YBa2Cu306+ x [0]- The peaks are observed in the magnetic response x"((f, uj) at Qi = (1±<5, 1)tt 
and Qi = (1, 1 ± 5)ir. It should be pointed out that some other theories such as the theory of 



stripes have natural explanations for the incommensurate peaks [BO, 100 1. 



3.5 Calculation of Vertex Correction 



Let us first assume that the quasiparticle residue is Z w 1. We will discuss the effect of the 



quasiparticle renormalization later. The first vertex correction is shown in Fig. |3.8| and can be 
written as 



«JA(3) = j J n(l - 2) G(l - 3) G(3 - 2) 
where numbers are representing the coordinates and LT is defined by 



(3.121) 



n 



< s{s{ >= - Yffi* - <rvV < s{s{ > 



3=1 



3=1 



a i \2<S\Si>-Y,<S{si> 

3=1 



-Ao X(l - 2) 



(3.122) 



where A D = (g/2)a l is the bare vertex. In the last step we assumed three equivalent modes of 
spin fluctuation, reflecting the rotational symmetry of the system. In momentum representation 
the first vertex correction can be calculated for general external momenta. Here we concentrate 
on the interaction between a Fermi surface electron and a spin fluctuation with q = Q. We 
have 



Ao 



d 2 q duj 
(2vr) 2 2^ 



X (Q' + q, w)G(k + Q' + q, u)G(k + Q" + q, w) 



(3.123) 
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k+Q , 




k +Q'+q CO 



Figure 3.8: First correction (5Ag to the bare vertex, for an incoming fcrmion with momentum k and 
energy (relative to the Fermi energy), interacting with a fluctuation of wave- vector Q and zero energy. 



where <5A^ = 5A(k,e = 0;q = Q,Q = 0) and Q" = Q + Q'; the sum over Q' takes care of 
Umklapp processes. In the reduced Brillouin zone Q" = (0,0). Note that the overall sign of 



the graph in Fig. 3.8 is negative, i.e., Eq. (3.123) is negative. 

We concentrate on SAr, with k being a point on the Fermi surface, since it gives an indication 
of the size of vertex corrections. This is because in general \5A(k, e; q, fi)| exceeds \6Ar\ (in fact 
it diverges along a surface in k-space, for given values of q, O and e); thus one cannot argue, even 
after integrating over one or more of its arguments, that |<SA(fc, e; q, 0)| will lead to corrections 
smaller than one would get from just using <5Ag (our argument here parallels Migdal's p6| ). 
More generally one finds that if 17 <C qvp,A , where A is the spin gap, then <5A^ is a good 
approximation to 5A(k, e; q, O). Writing ( |3.123| ) as 



Ao 



d 2 q duo x"(Q' + q,u) 



4 ' (2tt) 2 vr e 

Q' 



'k+Q'+q 



+ 



k+Q'+q 

i-/s 



'k+Q'+q 



-k+Q"+q 
fi 



k+Q"+q 



i-/s 



-g 



k+Q'+q 
2 XQ 



U) 



€ k+Q'+q +UJ 



z k+Q"+q 



U) 



e k+Q'+q +UJ , 



12vr 3 



LO s f fj, 



(3.124) 
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where g 2 = |<? 2 and is dimensionless; we define g to correspond directly with the coupling 



constant g used in Monthoux and Pines [35, |3(|. At zero temperature one has 

i , 2tv /■27T 

dq x i 

-2w J-2w " ei — €2 



/Z7T p2TT j 
dq x I dq y - — [Sgn(ei)Gi(|, ei) - Sgn(e 2 )Gi(f, e 2 ) - G 2 (q, ei) + G 2 (q, e 2 )] 
-2tt J-2tt ei — e 2 

(3.125) 



where ^ = g x a, etc., ei and e 2 are 

_ H+Q'+q ^ _ £ fc+Q"+q _ € k+q ,g 

1 A* ' 2 A* A* 

and G\ and G 2 are defined by 

G U9» e ) - 2 (e2 + X 2 ) ' G2 W' e )- g2 +X 2 (3.127) 

withX=[l + (i) 2 (g2+^)](u; s/ /|H). 

2r can be investigated both analytically and numerically. Here we calculate it for two 
different points in k-space, both on the Fermi surface (see Fig. |3.7| ); k\ makes a 45° angle with 
k x , and k^ is a "hot spot" wave vector. 

In order to obtain a value for SA^, we need values for the spin fluctuation energy oj s f, the 
correlation length £, the susceptibility \Q-> the coupling constant g, and the chemical potential 
\i (which is determined by the electron filling factor n). There are different values reported 
in the references [p3| , |32|. We have evaluated SA^ for k = k±,kh, in two ways, viz. (a) 
by assuming various published values for the different parameters, and (b) by making the 
simple assumption that one is very close to an AFM instability, and then, in the spirit of 
RPA, imposing the condition | (<?/2)x (0) (Q, 0)| « 1 where x (0) (£ 

u>) is the electron-hole bubble. 

Numerical calculation gives \x^(Qi 0)1 = 2.6(eV) _1 and thereby g = 0.67 eV, assuming n=0.75 
and the band structure in ( 3.112|) . Since uj s f/\/j,\ -C 1, this value of g should be a very good 
guess, within a naive RPA scheme. 

The results are summarized in Table 3.1. We use two different values for uj s f ; these are the 



two different values quoted by Monthoux and Pines et al. [35, [3j|. We also use two different 
values for the coupling constant g, quoted from Ref. [35] and [36]. We use values of £ = 2.5a, 
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9 (eV) 


uj s f (meV) 




SA R JA 




SAe/Ao 


MPI 


1.36 


7.7 


78.6 


-1.81 


105.6 


-2.43 


MPII 


0.64 


14 


49.6 


-0.46 


73.4 


-0.68 


"RPA" 


0.67 


7.7 


78.6 


-0.44 


105.6 


-0.59 


14 


49.6 


-0.50 


73.4 


-0.74 



Table 3.1: Calculated values of the vertex correction <5Ag for two different wave- vectors k\ and kh on the 
Fermi surface (columns 4 and 6 in the table). Ref. [§5l and [[36) give different values for g, and different 



values for w s f- From the values for these two models one calculates in equation (3.125), and thence 
5A^/A . The third model is the naive "RPA" model described in the text, for which g is determined; 
we have calculated and 8h.^/h. for two values of uj s j given in MPI and MPII respectively. 

Xq = 80 states/eV (from ]35|, |3(|) and n = 0.75 , appropriate to YBa2Cu3C>7 (again quoted 
from |35|, |36|). This value of n corresponds, with the band structure in ( |3.111| ), to a value of 
~ 1.46t = 0.365 eV. 

We see that even the values for the vertex correction calculated from the simple RPA model 
(b) are not small; as in the standard discussion of Migdal's theorem, the importance of vertex 



corrections appears in the ratio |<5A^|/A . If one takes values of g from the literature [35, 36], 
this ratio is quite unreasonably large (as large as 2.43 for the hot spots in the model used by Ref. 
[35]). Thus vertex corrections are clearly very important. The values we quote for \5A^\/A are 



considerably larger than previous estimates [23, 97, 101]. The reason for this difference with 
previous work can be tracked back to the factor I?, which is impossible to guess from purely 
dimensional arguments. In fact if we drop the factor 1^ from SA^, we get an order of magnitude 
estimate for 5A^ given by 

(3.128) 



oy»J^!^L\ « l 



A LJ 47T 3 |/x| 

which is broadly in agreement with previous estimates (see e.g. Millis [)97|]) ; x<2> u sf an d g 
must be redefined to conform with the parametrizations in this paper). 

In fact however is surprisingly large, and also shows a significant variation around the 
Fermi surface, with a maximum at the hot spots, and a minimum at intermediate wave-vector 
like ki . We should emphasize here that analytic calculations of 5A^ have to be approximated 
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rather carefully in order to give reasonable agreement with the numerical results in Table 3.1 
Approximations such as those of Hertz et al. p5|] (see also |101| ), which try to separate off 
a rapidly- varying (in cf-space) contribution from x"{q,uj), give quantitatively incorrect results 
(including a completely unphysical ln[(/c — kh)a] divergence as one approaches the hot spot). 

One might suppose that is large simply because of the band structure (i.e., because of van 
Hove singularities, or the hot spots). If this were true one could argue that the quasiparticle 
weight ought to be renormalized down near these singular points in the Brillouin zone, and 
that this would considerably reduce the vertex correction. In fact however we find this is not 
the case; this can be checked analytically by suppressing the regions immediately around the 
hot spots in the integral for Tr, or by simply redoing the numerical calculation for a slightly 
different band structure. We find that suppressing the hot spots entirely, reduces the vertex 
correction by a factor which is everywhere less than 2 (and which differs very little from unity 
when k is far from a hot spot). 

3.5.1 The Effect of Quasiparticle Renormalization Z 

Our calculation of the vertex correction was criticized by Chubukov, Monthoux and Morr 
(CMM) |5(J on the basis of inappropriate quasiparticle renormalization factor Z. As we men- 
tioned earlier the value we used was Z m 1 which is the value used in NAFL calculations of 
T c [31, 35, 36]. CMM's criticism was based on the fact that self consistent consideration of 
spin fluctuation damping gives spin fluctuation frequency uj s f in terms of other parameters. 



More precisely, one can write |)3| uj s f = c 2 w /2^ 2 j, where c sw is the spin wave velocity, £ is 
the correlation rate and 7 is the spin damping rate. In fact two points of views can be taken 
here. One can assume the spin environment as independent degrees of freedom. In that case 
u s f will be an independent parameter. On the other hand if one considers a one band model 
then assuming that the damping of the spin excitations is dominantly due to the spin-fermion 
interaction (and not other sources such as spin-spin exchange, impurity, etc.), the damping rate 
7 can be obtained from the imaginary part of the particle-hole bubble at transfer momentum 
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Q. The result of the calculation given in [50] is 



7T , , , V 2 



u sf = -\ shi(/>o|^=2 (3.129) 

4 g z Z 2 XQ 

where v is the Fermi velocity at the hot spots and cpQ is the angle between normals to the Fermi 



surface at the hot spots. CMM claim that if one uses the parameters reported in [36] and at the 
same time requires Z = 1, one finds uj s f ~ 1.06meV which is one order of magnitude smaller 
than Lo s j = 14meV used in Ref. [36] and our calculations. CMM then conclude that this is a 



sign of inconsistency in our calculations. Using the same parameters and solving (|3.129| ) for 
Z, with oj s f = 14meV one finds Z = 0.28. Substituting this value of Z in the calculations, 
the vertex correction will be reduced by a factor of Z 2 = 0.08. This way the magnitude of the 
correction will be one order of magnitude smaller than what we obtained, and therefore the 
vertex corrections will be small. 

We do not dispute the fact that all the calculations should be self-consistent. What we 
argue actually is that the self-consistency should be considered everywhere throughout the 
calculation. We also agree that the Hamiltonian ( |3.111| ) is valid for renormalized quasiparticles 
and therefore it is necessary to include the renormalization factor Z in all calculations. All 
the parameters we used in our calculation (including Z) were based on the reported values by 



Monthoux and Pines 36]. We emphasize again that the coupling constant g in Ref. [35 



is obtained by fitting the calculated T c to the measured value (~ 90K). If Z, as claimed in Ref. 



[50] is much smaller than one, then this value should be also considered in T c calculations and 
would affect the magnitude of g which gives a reasonable T c . 

In a strong coupling calculation, this means that all the fermionic lines have to be replaced 
by renormalized Green's functions. This multiplies the coupling constant g by a factor Z which 
exactly compensates the reduction of the vertex correction discussed in |3(]] . This is more clear 



by looking at ( 3.118| ). One realizes that a renormalized G R is going to reduce N(0) by a factor 
of Z. This factor multiplies g in the exponent in ( 3.11 7| ) . Keeping all the other parameters 
fixed the coupling constant should be increased to preserve the value of T c . The new coupling 
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constant is now given by 

9_ 
Z 



9 Z = ^ (3-130) 



where g is the (old) coupling obtained using Z = 1 (i.e. the values reported in Ref. [35] and 
||36|| ) . The same statement holds for weak coupling calculation ( |3.115| ) keeping in mind that 
A oc g 2 N(0) 2 . Substituting Z and g z into the calculation of the vertex correction we get 

(\ 5 H\\ 2 f \ 6A k\ \ 7 2 7 2 9 2 2 m 1q1 ^ 

K 3 — ► | — r 2 " I Z g z = Z ^ = g (3.131) 



The factor 1/Z in ( 3.130| ) therefore exactly cancels the effect of Z 2 in vertex correction cal- 



culations and our calculated value remains unchanged. Looking back at fl3.129j ) but now re- 
membering that gZ is actually what one obtains from fitting to T c , we find uj s j = 1.06meV 
even after the inclusion of the renormalization factor Z. The mismatch of this value with the 



reported value of Ref. [36] reflects the inconsistency of the theory rather than a mistake in our 
calculations. 

3.5.2 Sign of the Vertex Corrections 

The sign of the vertex correction is important because if it is positive, the vertex corrections 
tend to increase the effect of the spin fluctuations on the electron system resulting in a higher T c . 
On the other hand if it is negative, it decreases the spin-fermion coupling and reduces the T c . 
We saw that at the optimal doping, our calculation indicates a negative sign for the first vertex 



correction. On the other hand, Ref. [93, 50 1 report a positive sign at the optimal doping but 
negative sign near the antiferromagnetic instability. In this section we discuss what a change 
in the sign of the vertex correction between strong and weak doping would entail physically - 
essentially it implies a phase transition of some sort for which there is no experimental evidence 
that we know of. We give our argument based on a fact with which everyone agrees: near the 
antiferromagnetic instability the sign of the vertex correction has to be negative to suppress 



the effect of spin-fermion coupling 47] 



If the sign of the vertex corrections at the optimal doping is negative, as we found, then the 
vertex corrections tend to decrease the dressed coupling. Decreasing doping, increases the effect 
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of spin fluctuations and therefore increases the magnitude of the vertex corrections. This can 
be checked by Eq. ( p. 124 ). Decreasing doping would increase XQ (c* £ 2 ) an d thereby the vertex 



correction. Consequently, the dressed coupling would be decreased even further. This trend 
continues until the total vertex vanishes completely at the edge of the ordered phase. Thus the 
flow of the renormalized coupling when changing the doping is continuous. On the other hand 
if the vertex corrections have positive sign at the optimal doping, decreasing the doping would 
increase the dressed vertex. Continuation of this trend towards zero doping would give a very 
large dressed coupling which is unphysical. 



Interestingly it was speculated by Chubukov et al [51, 93 1 that perhaps a weak phase tran- 
sition (a Lifshitz transition ) might exist in the intermediate doping regime. The transition 
would change the big Fermi surface, compatible with Luttinger's theorem, to small hole pock- 
ets centered at (7r/2a, 7r/2a) and symmetry related points. If this were really true, and moreover 
followed from the Hubbard model, it might save the consistency of the NAFL approach, in that 
it would make it possible for a sign change of the vertex corrections, at some critical g = g c . 
To our knowledge, no experimental justification exists yet ]52[ and the evolution of the Fermi 
surface is not widely believed. Without this evolution, the positive sign of the vertex correc- 
tion at optimal doping is not physical. At the end we point out that the suppression of the 
spin-fermion coupling due to the vertex corrections, also gives a natural way to explain why T c 
reduces at low dopings in spite of the increase in the spin susceptibility (9^] (although other 
reasons can also be thought of, e.g. reduction of charged carriers, etc.). 

3.6 Concluding Remarks 

We re-emphasize here that these results do not necessarily invalidate the internal consistency 
of the Fermi liquid starting point, in NAFL theory. However they do show that the theory 
cannot be trusted quantitatively, at least in the usual RPA form. As is well known the RPA is 
not a "conserving approximation" , and for spin fluctuation theories this makes it unreliable (cf. 
ref.(6l|], especially section 3). It is useful to compare the case of nearly ferromagnetic 3 He liquid, 
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where vertex corrections are also quite large, and where use of the paramagnon model yields 
values for m* /m which are off by a large factor J102j ] . Thus if we use melting curve Landau 
parameters, Zq ~ 0.75 and F\ ~ 15, we infer a value for the Stoner factor S ~ 24 which yields 
m* /m = | InS ~ 15 , in the paramagnon model. This is roughly 2.5 times the correct value of 
~ 6 (note that the first vertex correction is 5 A/ 1 ~ In S ~ 3 in this model), and no amount of 
self-consistent summing of diagrams can cure this numerical problem. 

Similar problems can clearly occur in the present NAFL model. We believe this is the main 
reason for the difficulty one encounters in the MP models, in determining a value for g that (a) 
gives the correct superconducting T c , and (b) is consistent with the observed spin susceptibility. 

To check the structure at higher order, we have also estimated the contributions from the 
graphs containing 2 spin fluctuation lines (there are actually 7 distinct graphs at this level) , and 
found that some of them are also large for the values of g used above. Thus, just as for the case 
of nearly ferromagnetic 3 He, we see no reason to believe, for the values of the parameters given 
in the table, that performing infinite graphical sums will lead to results which are numerically 
more reliable, even if they do converge to some smaller renormalized vertex-there will always 
be other diagrams with large values, which will in general give uncontrolled contributions. 

It is interesting to compare our results with some other investigations. As we mentioned 



briefly, in the weak-coupling limit, Chubukov [49, 93 1 has calculated the leading vertex correc- 
tions to g, concentrating on the gapless case; in the case where there is a gap, he finds that the 
renormalized coupling is large (for a large Fermi surface considered in our investigation). On 



the other hand Schrieffer [47] has argued that a correct formulation of the theory, even in the 
weak-coupling limit, must take account of the short-range local antiferromagnetic order even 
in the normal state-if this done, he finds that a weak-coupling calculation shows very strong 
suppression of the vertex when one is close to the antiferromagnetic transition. This theory 



seems rather interesting-note that a related calculation by Vilk and Tremblay [64, |103[| finds 
that the existence of such a short-range antiferromagnetic order in 2-d will cause a breakdown 
of the Fermi-liquid starting point itself when the correlation length becomes larger than the 
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single particle de Broglie wavelength! Thus the question of what is the correct theory itself is 
rather confused, even in the weak-coupling regime. It is certainly not clear how any of these 
arguments will work in the regime discussed in this paper, when g is not small enough to control 
the magnitude of the vertex corrections. 

It is of course crucial that these higher-order corrections also be included in any version 
of this theory that tries to reconcile different experiments - as emphasized by Pines |32|], the 
justification of the theory stands or falls on its ability to do this quantitatively. It is possible 
that such a program might succeed if one can show that the actual parameters g, cj s f, and XQ 
are such that |<5Ar|/A is considerably less than one (i.e., if one is genuinely in the weak-coupling 
regime). This would also be true of versions of the theory in which u) s f depends on g, whilst 
the spin gap becomes an independent parameter |^|; or of the theory of Schrieffer cited above 
pTfl . On the other hand if |<5A^| / A > O(l), we see no hope that such a scheme could succeed 
quantitatively (in, e.g., the calculation of T c ), since the vertex corrections become large. 
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Chapter 4 



Introduction 



Although at a microscopic level, theories of high T c superconductors, especially for normal 
state properties, are not well established, superconducting properties of cuprates are better 
understood. Almost everybody now agrees that the symmetry of the order parameter is d x 2_ y 2- 
wave, and the existence of nodes in the superconducting gap is well established |33|, [Mj. A 
highly anisotropic d^.^-wave order parameter, with four nodes in the superconducting gap 
is a unique property of cuprates ans also some organic superconductors, with no analog in 
conventional superconductors. The effect of this anomalous symmetry on different observable 
quantities has been a subject of investigation in recent years. 

An early theoretical investigation of the weak-field response of a d x 2_ y 2 superconductor 
by Yip and Sauls [104] predicted a direction dependent non- linear Meissner effect, associated 
with the quasiclassical shift of the excitation spectrum due to the superfiow created by the 
screening currents. Maeda et al. |105|1 reported experimental evidence for such an effect in 



Bi2Sr2CaCu20 y , but subsequent experiments [106] failed to confirm their findings and the sit- 



uation remains controversial. A similar effect was also studied independently by Volovik [107]. 
In the mixed state, he predicted a contribution to the residual density of states (DOS) propor- 
tional to the inter-vortex distance ~ \f~H. Such a contribution was identified in the specific 



heat measurements on YBa2Cu30y_5 (YBCO) by Moler et aZ.[108] but later this interpreta- 



tion was disputed by Ramirez [108] who found evidence for a similar effect in a conventional 
superconductor VsSi and by others[llC]. Kosztin and Leggett ]|111| predicted that the nonlocal 
response at very low temperatures will lead to a T 2 dependence of the penetration depth in 
clean samples in contrast to the linear T-dependence obtained from the local theory for d-wave 
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materials. 

Another important characteristic of high T c compounds is their short coherence length £o> 
responsible for strong fluctuation effects near T c [|112| . The large ratio k = Ao/Co 3> 1 (Ao is the 
magnetic penetration depth) also makes them extremely type-II superconductors. Existence of 
a mixed state, characterized by a regular array of magnetic flux lines penetrating the material, is 
perhaps one of the most striking properties of type-II superconductors. The original pioneering 



work of Abrikosov [113], based on the solution of Ginsburg-Landau (GL) equations [114] near 
the upper critical field H C 2, predicted a triangular flux lattice. This prediction was subsequently 
verified by low field magnetic decoration experiments on a variety of conventional supercon- 
ductors. In stronger fields neutron scattering experiments revealed significant deviations from 
perfect triangular lattices [115] which where attributed to anisotropies in the electronic band 
structure and other effects. These were modeled by GL theories containing additional higher 
order derivative terms reflecting the material anisotropies |115f |. 

Based on the experience with conventional superconductors one would expect even richer 
behavior of flux lattices in the copper-oxide superconductors. In high-T c cuprates much of the 
experimental and theoretical effort has been focused on the sizable region of the phase diagram 
just below T C (H) in which the vortex lattice properties are completely dominated by thermal 
fluctuations El , While understanding the physics of this fluctuation dominated regime poses 
a very difficult statistical mechanics problem, investigation of the equilibrium vortex lattice 
structures at low temperatures may provide clues about the microscopic mechanism in these 
materials. 



Neutron scattering [117] and STM [118] experiments on YBa2Cu30"7_,5 compound, revealed 
vortex lattices with centered rectangular symmetry and different orientations with respect to 
the ionic lattice. These have been modeled by phenomenological GL theories appropriate for 



anisotropic superconductors, containing additional quartic derivative terms [119] or a mixed gra- 



dient coupling to an order parameter with different symmetry [12C, 121, 122, 123]. These works 
found structures in qualitative agreement with experiment, but their inherent shortcoming is 
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a large number of unknown phenomenological parameters and subsequent lack of predictive 
power. Also, the GL theory is only solvable for a vortex lattice near H C 2, which is experimen- 
tally inaccessible in cuprates away from T c . The observed behavior of the vortex lattice may 
also be understood by incorporating penetration depth anisotropy and twin-boundary pinning 



without involving any effects associated with mixing symmetries or gap anisotropy [124]. 



Muon-spin-rotation (/iSR) experiments |125j , |126| , |127 , 128 ], on the other hand, show an 



unusual magnetic field dependence in their line-shapes for the magnetic field distribution. This 
has been attributed to a field dependent penetration depth - which is expected in the Meissner 



state because of quasiparticle accumulation at gap nodes |104| . It was also modeled using an 
approach based on the Bogoliubov- de Gennes (BdG) equations in a square lattice tight-binding 
model (129 1 . 



At intermediate fields H c \ < H <C H C 2, properties of the flux lattice are determined primar- 
ily by the superfluid response of the condensate, i.e. by the relation between the supercurrent j 
and the superfluid velocity v s . In conventional isotropic strongly type-II superconductors this 
relation is to a good approximation that of simple proportionality, 

j = -ep s v s , (4.132) 

where p s is a superfluid density. More generally, however, this relation can be both nonlocal 
and nonlinear. The concept of nonlocal response dates back to the ideas of Pippard |130| 1 and 
is related to the fact that the current response must be averaged over the finite size of the 
Cooper pair given by the coherence length £o- I n strongly type-II materials the magnetic field 
varies on a length scale given by the London penetration depth Ao, which is much larger than 
£o and therefore nonlocality is typically unimportant unless there exist strong anisotropics in 
the electronic band structure [ |131| . Nonlinear corrections arise from the change of quasiparticle 
population due to superflow which, to leading order, modifies the excitation spectrum by a 
quasiclassical shift |13C|1 

S k = E k + v r v s , (4.133) 
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where Ek = y e| + is the BCS energy. Again, in clean, fully gapped conventional supercon- 
ductors this effect is typically negligible except when the current approaches the pair breaking 
value. In the mixed state this happens only in the close vicinity of the vortex cores which oc- 
cupy a small fraction of the total sample volume at fields well below H C 2- The situation changes 
dramatically when the order parameter has nodes, such as in d x 2_ y 2 superconductors. Nonlo- 
cal corrections to ( [4.132 ) become important for the response of electrons with momenta on the 
Fermi surface close to the gap nodes even for strongly type-II materials. This can be understood 
by realizing that the coherence length, being inversely proportional to the gap |130(| , becomes 
very large close to the node and formally diverges at the nodal point. Thus, quite generally, 
there exists a locus of points on the Fermi surface where £ ^> Aq and the response becomes 



highly nonlocal. This effect was first discussed by Kosztin and Leggett [111] for the Meissner 



state and by us [ 132 1 in the mixed state. Similarly, the nonlinear corrections become important 
in a d-wave superconductor. Eq. ( 4.133| ) indicates that finite areas of gapless excitations appear 
near the node for arbitrarily small v s . 

One of the simplest theories to study the magnetic behavior of the type-II superconductors is 
the London model | |133j ]. In London theory the Free energy is written only in terms of magnetic 
field and is therefore very easy for calculation. Unlike in GL theory, the order parameter 
does not appear explicitly in London free energy. The effect of the symmetry of the order 
parameter or other Fermi surface anisotropies, is therefore not contained in the London model. 
Simple generalization of the model however can incorporate these effects. The generalized model 
contains anisotropic higher order and higher derivative terms reflecting the symmetries of the 
order parameter and the anisotropic Fermi surface. 

In the next chapter, we first briefly introduce Ginsburg-Landau (GL) theory, and London 
theory emphasizing the solutions for a vortex lattice. Starting from a GL model, we then derive 
the leading fourfold anisotropic corrections to the London equation making the usual assumption 
that the free energy was an analytic functional of the order parameter and field. The number 
of new parameters in this model is far smaller than in the GL approach (a reasonable model 
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contains only one new parameter which controls the strength of the symmetry breaking term) 
and numerical simulations are considerably easier. This provides a useful tool to study vortex 
lattice structure, pinning by twin boundaries and the magnetic field distribution measured in 
fiSR experiments. The model is suitable to study the intermediate field region H c \ <C H <C H C 2 
which is experimentally most relevant but traditionally difficult to handle within the GL theory. 
Furthermore, this approach can be extended to T = where GL theory breaks down and the 
supercurrent becomes singular. With increasing magnetic field this model predicts a transition 
from triangular to centered rectangular and eventually a square vortex lattice. While no direct 
experimental evidence exists in cuprates at present to confirm such a prediction, a similar 



transition has been recently observed in a boro-carbide material ErN^E^C [134] and has been 
described by a similar London model 131 1. 

In the last chapter of this thesis, we again derive a generalized London free energy, but 
this time completely from a weak coupling microscopic model, including the nonlinear and 
nonlocal effects mentioned above. We show that the dominant effect that determines the vortex 
lattice geometry and the effective penetration depth as defined in /iSR experiments is the 
nonlocal corrections, while the nonlinear corrections play a secondary role at low T. At high 
temperatures we obtain a nonlocal correction similar to the one suggested in chapter 5. At low 
temperatures however, we find a novel singular behavior directly related to the nodal structure of 
the gap which completely changes the form of the London equation. This singular behavior has 
profound implications for the structure of the vortex lattice which, as a function of decreasing 
temperature, undergoes a series of sharp structural transitions and attains a universal limit at 
T = 0. Our theory is now completely parameter free. The London free energy at low T is 
non-analytic and its long wavelength part is fully determined by the nodal structure of the gap 
function. Such behavior is caused by the low-lying quasiparticle excitations within the nodes 
and thus could never occur in conventional superconductors with anisotropic band structures. 

Finally, in the last section of chapter 6, we talk about the experimental justification of our 
low temperature nonlinear nonlocal theory. We especially emphasize //SR experiments. Recent 
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fiSR data at high fields have justified our prediction for the Field dependence of the penetration 
depth in vortex state. 



Chapter 5 



Phenomenological Model 



5.1 Ginsburg-Landau Theory 



In 1950 Ginsburg and Landau [114] proposed a phenomenological theory to describe the be- 
havior of superconductors near their transition temperature. At a mean field level, the phase 
transition of superconductors is second order (gauge field fluctuations however are believed to 
change the order of the transition to a weak first order |135fl ). A free energy density that can 
describe such a transition was proposed to be 

/ = -^K-i/iV - -A)^\ 2 + aM 2 + ^| 4 + (5.134) 
2m* c 2 57T 

where e* (= 2e) and m* (= 2m) are the charge and effective mass of a Cooper pair, H and 
A are magnetic field and vector potential respectively Superfluid density n s is related to the 
order parameter by l^l 2 = n s /n, where n is the total density of electrons. Ginsburg-Landau 
(GL) equations can be obtained by minimizing / with respect to ^ and A 

a* + pm 2 ^ + — \(-ihV - —A) 2 ^ = (5.135) 
2m c 

J = -^-V x B = - — (f* W - fW*) - —y*^A (5.136) 

4tt 2im* m*c 

where j is the current density. In the absence of magnetic field (i.e. H = A = 0) and for 



uniform ^, the minimum of / occurs at \^\ =0 when a > 0, but at |^| = y/—ot//3 when a < 0. 
A continuous transition therefore happens as a goes from positive to negative values. Taking 
a oc (T — T c ), one can study the critical behavior of the superconductor and calculate all the 
critical exponents near T c . GL theory has been studied extensively in the literature in both 
Meissner and vortex states of superconductors. It is not our attempt in this thesis to discuss the 
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solutions of GL equations for different problems. Instead, we shall focus on deriving a London 
theory from GL theory, especially in the vortex state. Interested readers should refer to Ref. 
fll36j , |130| 1 or any other standard text books in the subject, for more detailed information. 

5.2 London Theory in Vortex State 

Although we will use GL equations to derive the London equation, London theory was proposed 



long before GL, by F. and H. London in 1935 [133]. The basic assumption in their theory is a 



local proportionality relation between supercurrent j and vector potential A 

-A. (5.137) 



n,e 2 r 



mc 



This equation is usually taken as the definition of superfluid density n s . Eq. (5.137) is not a 
gauge invariant equation. However, to satisfy current conservation, one has to fix the gauge to 
the London gauge 

V-A = 0. (5.138) 
This way ( |5.137 ) and all the results following from it are gauge independent and therefore 



physical. Eq. (5.137) plays a very important role in our theory for the vortex lattice. General- 
ization of this equation will be used in this chapter and also in the next chapter to develop a 
generalized London theory to describe the vortex lattice properties of d-wave superconductors. 

Taking the curl of the both sides of ( |5.137| ) and using the Maxwell's equation V x B = 
(47r/c)j we get 

A 2 VxVx6 + B = (5.139) 
which is the famous London equation. A is called the London penetration depth and is given by 



mc 2 



< 5 - 140 > 

A free energy density that gives rise to ( |5.139| ) via minimization with respect to the magnetic 
field B, is called the London free energy density 

f L = ^(B 2 + X 2 \VxB\ 2 ) (5.141) 
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We now obtain London equations ( |5. 139 ) from GL equations ( |5.135| ) and ( |5,136j ). 



5.2.1 Derivation of the London Equation from GL Theory 

Let us begin our derivation of the London equation by writing the order parameter as \l/ = ipe 1 ^ , 
with ift and eft being real functions of position x. Substituting ^ into ( |5.136 ) we get 



\ 2 V x B + A = p-V<f) (5.142) 



where A = \J m*c 2 /Ane* 2 ift 2 is the London penetration depth and $o = (hc/e*) is the flux 
quantum. In type II superconductors ift varies over a length scale £ much smaller than the 
penetration depth A. £ is usually called the coherence length and is of the order of the extent of 
a cooper pair. In the Meissner state therefore, tp (and consequently A) stays constant in the bulk 
of the superconductor away from the boundaries. The London equation ( 5.139| ) can therefore 



be derived easily by taking the curl of Eq. (5.142) and keeping in mind that V x Vcft = 



In the vortex state on the other hand, ift is not a constant and actually vanishes at the center 
of each vortex. Taking this into account gives rise to the effect of the vortex core which is the 
subject of the next subsection. Furthermore, in the derivation of ( |5.139| ), we used V x X7<ft = 
which is not true in the vortex state. The reason is that around a vortex, (ft has a topological 
winding number equal to 1. In other words, (ft increases by 2tt in every turn around a vortex 
core. Integrating (V x V0) • dS over a small region around the vortex center we get 

J (V x Vcft) ■ dS = j>V(ft- dl = <ft f -<Pi = 2vrn (5.143) 

where n is the topological winding number; which is one in our case. This equation is satisfied 
for every infinitesimal region around the center of the vortex ro, only if V x = 27rj( 2 ) (r — tq) . 
The London equation in the presence of a vortex will therefore become (assuming constant tft 
again) 

A 2 VxVxB + J] = f/)(f-f ) (5.144) 
To find a solution for a single vortex, it is easy to take the magnetic field in z-direction and 
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write B = B(r)z. Eq. ( |5.144| ) then becomes 

" — T-C'^r) + B = <M (2) (r- r ) (5.145) 
r dr dr 

which is actually a Bessel's differential equation. The solution to this equation is 

where Kq is the zero-order modified Bessel function. Important to notice here is that Kq 
diverges logarithmically as r — > 0; which is apparently non-physical. The reason behind this 
non-physical behavior lies in the fact that we did not incorporate the effect of the vortex core 
properly. In other words, the order parameter tjj always vanishes at the center of the vortex. 
Therefore, it is not a good approximation to consider it to be homogeneous everywhere in the 
lattice. We shall take into account this effect in more detail in the next subsection. 

5.2.2 The Effect of the Vortex Core 

Let us start from eq. (5.136) again by writing it as 



j{r) = e*n s (r)v s (r) (5.147) 
where n s = if) 2 is the superfluid density and v s is the superfluid velocity defined by 

At the center of a vortex, n s vanishes. We can therefore write 

n s (r) = norj(r) (5.149) 
where rj(f) is a function with the property 

riir) = at r = 

(5.150) 

rj(r) — ► 1 for r/£ 3> 1 

and £ is the coherence length of the superconductor. The exact form of rj should come from a 
microscopic theory or a complete solution of GL equation. Here, we do not derive rj. Instead, 
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we find the source function for a few assumed forms for 77. Before that, let us take the curl of 
both sides of ( p.!47| ), keeping in mind that 



V x (r)v s 



2-kTi 



m 



-7?(<5 (2) (r)--^B) + V77.. . s 



x v s = — —^B + vrj x v s 



and also from (5.147) 



1 



J 



The generalized London equation therefore becomes 



V x B 



A 2 V x V x B + B = p(r) 



with a source term p{r) defined by 



p{f) = (1 - n)B + A — - x V x B 
V 



(5.151) 



(5.152) 



(5.153) 



(5.154) 



Equations ( |5TT53|) and (|l5|) can be viewed as a set of self-consistent equations to be solved 
for B and p. However, p(r) can be simplified significantly in extremely type II superconductors 
which have large k = A/^. First notice that p is non-zero only in a region of size £ around the 
center of the vortex. In this region, both B and r] scale with £. Therefore the first term in 



the left hand side of (5.153) and the second term on the right hand side of (5.154) are both 
0{\ 2 /i 2 )B and therefore are the dominant terms in ( |5.153| ), Ignoring the subdominant terms, 
(|5.153[) can be written as 



V x u = — x u 



V x 







(5.155) 



7] \7] 

where u = V x B. If we assume rotational symmetry and let B = B(r)z and therefore u = u(r 
where 6 is the azimuthal angle, then ( |5.155| ) leads to 



Id / u 
r dr \ rj 







u = C- 



(5.156) 



C is a constant that can be found, using the fact that J p ■ dS = &o, to be C = (<3?o/27rA 2 ). 
Substituting all into ( |5.154 ), the source function now can be written as p(r) = p(r)z, with 

$0 drj 



p(r) 



2nr dr 



(5.157) 
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Now let's calculate p for different choices of rj. 



Example 1: Clem's Ansatz 

Let's assume that 

rj(r) 



r 

R 2 



with 



R 2 = r 2 + £ 2 



(5.158) 



where £ v is some parameter of the order of the coherence length £. Substituting into ( |5,157j ), 
we find the source term to be 



p{r) 



vr R A 



(5.159) 



In order to calculate the magnetic field, it is easier to find the source term in terms of the 
magnetic field, using ( |5 . 1 54 ) . Ignoring the first term we get 

2H dB{r) 



p{r) 



R 3 dR ' 



(5.160) 



Substituting back into (|5.153|) , with a little algebra, we get 

R dR V dR) 
Comparing ( |5,161 ) with ( 5.145|) , the solution would be 



(5.161) 



B(r) 



2vrA 2 



(5.162) 



This result was first obtained by Clem [137], directly by solving GL equation. It was later 



shown by Brandt [138] that the magnetic field obtained this way in a vortex lattice is consistent 
with the exact solution of the GL equations. In Clem's approach, £ v is a variational param- 
eter to be determined by minimizing the free energy. It turns out that in extremely type II 
superconductors = v2£- As is clear from ( |5.162 ), B(r) is not divergent as r — > 0. In fact, 
B(0) = (<Do/2^A 2 )A^/A). 

Another way to find the magnetic field is by Fourier transforming ( [5.1 53j ) to get 

d 2 k F{k)e^ 



\ 2 k 2 B{k) + B{k) = F{k) or B{r) 



(2tt) 2 1 + X 2 k 2 



(5.163) 
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where 



F(k) = I d 2 re~ ik - r p(r) 



(5.164) 



is the Fourier transform of the source function which sets an ultraviolet cutoff in momentum 
integrations. In the present case, it is not difficult to show that 



F(k) = *o&fctfi(&A:) 



(5.165) 



where K\ is the first order modified Bessel function. Substituting ( |5.165| ) into ([5.163 ), we get 



back to (|5l62D . 



Example 2: Gaussian Source Function 

Now let us assume that 

r,(r) = 1 - e" r2 /2« 2 



The resulting source function will then be 



$o drj = $ c -r 2 /2g 2 



2nr dr 2tt^ 2 

The corresponding cutoff function will also have a Gaussian form 

F(k) = cD oe -« 2fc2 / 2 



(5.166) 



(5.167) 



(5.168) 



This form of momentum cutoff was first proposed by Brandt [13£]. This is actually the form 
of cutoff we are going to use in our calculations throughout the rest of this thesis. We will 
talk about the sensitivity of our solutions to this particular choice of cutoff function, when 
appropriate later. 



5.2.3 Solution for a Vortex Lattice 



In a vortex lattice, taking the magnetic field in z-direction {B = Bz), equation ( |5. 153 ) becomes 



X 2 V 2 B + B = J2p( 



r — n 



(5.169) 
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where fj are lattice vectors indicating the position of vortices. Taking the Fourier transform of 
both sides of fl5.16Sp , we get 



1 + \ 2 k 2 



B(r) 



„ik-f 



F(k) 



+ \ 2 k 2 



(5.170) 



where k are now reciprocal lattice wave vectors, B is the average magnetic field given by 
B = &o/£l, with f2 being the area of a unit cell. Unlike GL, the London approach does not 
automatically determine the form of the lattice for us. Instead, one has to calculate the free 
energy for different configurations of the lattice and find the one that minimizes the free energy. 
Although this sounds a formidable job, symmetries of the free energy help us in our initial 
guess. In other words, the lattice that satisfies the symmetries of the free energy the most, is 
the most likely lattice for the vortices to exist on. 
















Figure 5.9: Magnetic field distribution in a triangular lattice, using the ordinary London equa- 
tion with Gaussian cut off. 



+ X 2 k 2 



Using ( [5.141| ), the free energy per unit cell is given by 

f L = [ d*r[B 2 + (VBf\ = £(1 + X^)Bl = B*Y: 
Juc ~ ~ 1 

k k 

Substituting the Gaussian cutoff function ( |5.168| ), we find 

p ik-f p -^k 2 /2 



(5.171) 



(5.172) 
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e< 2k 



k 



These sums can be performed numerically quite easily. The lattice configuration that minimizes 
the free energy (|5.173 ) is a triangular lattice. The contour plot of the magnetic field in this 



lattice is shown in Fig |5.9| . 

An important thing to notice is that the free energy ( |5.173| ) is isotropic and therefore does 
not determine the orientation of the lattice. In other words, there is an infinite degeneracy in 
the orientation of the lattice. This degeneracy will be removed as soon as symmetry breaking 
terms are added to the London equation. We will discuss this issue in more detail in the next 
two sections. 

5.3 Mixing of d-wave with s-wave 

So far we have established a London theory derived from GL theory, but we didn't mention 
anything about the symmetry of the order parameter. Since an s-wave order parameter is 
isotropic, whatever we discussed in the previous sections is expected to hold for an s-wave 
superconductor. A d-wave order parameter on the other hand, is not isotropic. Instead, it 
has fourfold anisotropy, and consequently, it is natural to expect this anisotropy to influence 



the structure of the vortex lattice. Surprisingly, the GL equations, as written in ( 5.135 ) and 



(5.136), have exactly the same form for a d-wave order parameter as s-wave. All the statements 



in the previous sections should therefore hold for a d-wave superconductor as well. A way 
out of this dilemma, is to notice the fact that the GL free energy written in ( |5.134 ), is an 



approximation and contains only the lowest order contributions. In other words, we have 
ignored all higher order derivative terms. We have also overlooked the possibility of existence 
of a second order parameter that in general can couple to the main (d-wave) order parameter. 
It is conceivable that including these effects will introduce some fourfold anisotropic terms to 
the GL equations, reflecting the d-wave nature of the superconductor. In this section, we shall 
first consider the second possibility, namely including an s-wave order parameter coupled to the 
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original d-wave order parameter. We shall see that the effect of this inclusion is to introduce 
higher order derivative terms in the effective GL free energy written only in terms of the d- 
wave order parameter. We shall then show how this additional high-derivative term will change 
the London free energy. This way we actually cover both possibilities mentioned above. The 
material discussed in this section has been published in Ref. |140| ] . 

A GL free energy density to describe mixing of d-wave and s-wave order parameters was 
first proposed by Joynt |141| 1 as 



/ = a s \s\ 2 + Ud\d\ 2 + "f s \ns\ 2 + ld \nd\ 2 + U + H 2 /8ir 

+ lv [(U y s)*(U y d) - (U x s)*(U x d) + c.c.]. (5.174) 

where II = —iSJ — e* A/fic and contains the quartic terms, d and s represent d-wave and 
s-wave order parameters respectively. In the bulk of the superconductor and in the absence of 
magnetic field, we want only a d-wave order parameter to exist. Therefore below T c , we expect 
to have a d < 0, but a s > 0. In finite field (H > H c \) on the other hand, a small s-component 
with a highly anisotropic spatial distribution is nucleated in the vicinity of a vortex [120, 121 



This free energy has been extensively studied in the literature 120 , 121 , 122 , 123 , 142]. It is 



well known that it gives rise to non-triangular equilibrium lattice structures [120, 121 1 . Our 
objective here, is not to discuss the solution to the GL equations resulting from this free energy, 
but to focus on deriving the corrections to the London equation that result from this free energy. 



Our strategy will be to simplify the free energy (5.174) by integrating out the s-component in 
favor of higher order derivative terms in d. In this process, some short length-scale information 
on the order parameter is lost but the magnetic field distribution is described accurately. 

Using its Euler-Lagrange equation, s can be expressed to the leading order in (1 — T/T c ) as 

s = (7„/a a )(Il£ - U 2 y )d. (5.175) 

Substituting this into / gives the leading derivative terms in d of the form: 

/ = ld [\fld\ 2 - ( 7 2 / ld a s )\(Ul - U 2 )d\ 2 ] + ... (5.176) 
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Various additional corrections to the free energy are obtained from integrating out s more 
accurately, taking into account the 7 s |IIs| 2 term and quartic terms. However these all involve 
higher powers of II or other terms that will not concern us. The coefficient of the second 
term has dimensions of (length) 2 . For future convenience, we will write it in the form e£ 2 /3, 
where e = 3(a ( i7^/Q ; s7j) is a dimensionless parameter which controls the strength of the s-d 



coupling and £ = vTd/l^dl is the GL coherence length. We henceforth assume e <C 1. As 
we will discuss in the next section, neutron scattering and STM experiments probably support 
this assumption. As mentioned before, a term of the form |(II 2 — II 2 )<i| 2 could arise without 
invoking s-d mixing. It can result from a systematic derivation of higher order terms in the GL 
free energy starting with a BCS-like model and taking into account the square symmetry of the 



Fermi surface p|, flm pH, 145 1 • 



The free energy ( |5.176 ) is not bounded below, exhibiting runaway behavior for rapidly 
varying d-fields. This is in fact cured by keeping additional higher derivative terms that also 
arise from integrating out s. In fact, the approximation of Eq. ( |5,176| ) will be sufficient for our 
purposes, yielding a local minimum which we expect would become a global minimum upon 
including the additional terms. 

In the case of extremely type II superconductors - which we are considering here - the 
penetration depth A ^> £. We may therefore assume that \d(r)\ ~ do, the zero field equilibrium 
value, almost everywhere in the vortex lattice, except within a distance of 0(£) of the cores. 
This gives the London free energy, 

f L = (l/8n)(Bf + ld dl{v 2 - (ee/3)[(v 2 x - + {d y v y - d x v x ) 2 ]}, (5.177) 

written in terms of the superfluid velocity (we set fi = m* = 1), 

v = X79 - (e*/c)A, (5.178) 

where 9 is the phase of d. 

The corresponding London equation, obtained by varying fi with respect to A, is: 
c -> /2e*\ 2 

—V x B = [j^j lddl{v- -ef[(yv y - xv x )(v 2 y - t; 2 ) - (yd y - xd x ){d y v y - d x v x )}}. (5.179) 
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For many purposes it is very convenient to express v in terms of B and its derivatives, and then 
substitute this expression for v back into /l, giving an explicit expression for /l as a functional 
of B only. For e = this gives 

tfW = V x B/B , (5.180) 
where Bq = <&o/2ttAq is a characteristic field of order H&, and 

/E = (l/87r)[S 2 +Ag(VxS) 2 ]. (5.181) 



which is actually the ordinary London free energy (i.e. Eq. ( 5.141 )). Here the penetration 



depth, for e = is A 2 = 87r^d{e* do / he) 2 . It is presumably not possible to solve Eq. ( 5.17S| ) 



in closed form for v as a function of B for t / 0. However, this can be done readily in a 
perturbative expansion in e. The first order correction is: 

= (2ef/3){(yv^ - xv®)[(v®) 2 - (v x 0) ) 2 ] - (yd y - xd x )(d y v^ - d x v^)} (5.182) 

with given by Eq. (|5l80| ). The London free energy density, up to 0(e) is then: 



h = f°L+ € -^[4(d x d y B) 2 + ((d x B) 2 - {d y B) 2 ) 2 /B 2 ]. (5.183) 

Note that we could have arrived at a similar conclusion by simply writing down all terms allowed 
by symmetry in fi, expanding in number of derivatives and powers of B. Square anisotropy is 
first possible in the fourth derivative terms. In principle, we should also include all isotropic 
terms to order B 4 and V 4 . However, assuming that these have small coefficients, they will not 
be important. As we shall see in the next chapter, similar results can also be obtained from 
nonlocal effects due to the divergence of the coherence length along the node directions [132] 



and also from considering generation of quasi-particles near gap nodes | |104j , |146| ], in a range of 
temperature and field where the supercurrent can be Taylor expanded in the superfluid velocity. 
More generally, the quadratic and quartic terms in (5.183) have independent coefficients. 



The corresponding London equation is obtained by varying fi with respect to B(r). For B 
along the z direction, one obtains 



[1 - X 2 V 2 + Ae\lf{d x d y f)B - eQ[B] = 0, (5.184) 
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where 



Q[B] = 2A^ 2 VKd* " d y) B + d xBd x - d y Bd y ][(d x B) 2 - (d y B) 



(5.185) 



is the non-linear term arising from the last term in Eq. (5.183). 

Our numerical calculation (next section) shows that, the effect of the non- linear term is 
negligible (Contrary to naive expectation, it doesn't become more important with increasing 
applied field because the field becomes nearly constant in the vortex lattice when the applied 
field is large.) Thus to an excellent approximation one may neglect Q[B] in the London equation 



( 5.184 ). To get a feeling for the effect of the extra nonlocal (four derivative) term, consider 
a semi-infinite superconductor in its Meissner state with an interface parallel to the yz plane. 
The magnetic field in this case depends only on x; all the y-derivatives and thereby the extra 
term vanish. The penetration depth will then remain unchanged and equal to Ao- On the 
other hand, if the interface is 45° rotated with respect to x and y axis (o and b directions), the 
magnetic field will be a function of (x + y). Substituting B = Boe~^ x+y ^^ x into ( 5.184 ), we 
get (neglecting the nonlinear part) 

1/2 



A = Ar 



1 



1 



(5.186) 



Here we keep only the larger solution because it is the one that determines the decay rate of 
the field and therefore can be identified as the penetration depth. The penetration depth is 
now longer along the crystal axes as is clear from (5.186). 

To determine vortex lattice structure, we insert source terms J2j P(^~^j) a t the vortex core 



positions, fj, on the right hand side of Eq. (5.184). The source term we use is the Gaussian 
source term introduced in eq. ( 5.167| ). Taking Fourier transform, the magnetic field may be 
written explicitly as: 

8| '>^IWTi«W' (5 - 187) 

k 

Here the sum is over all wave-vectors in the reciprocal lattice and B is the average field. 
The lattice constant is determined by the condition that BQ, = <3?o, where f2 is the area of 
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the unit cell. The lattice symmetry is then determined by minimizing the Gibbs free energy 
Gl = J~L — HB/Air, where 

^ L = B Vi + WT^mkM 2 (5 - 188) 

k 

5.4 Numerical Results and Experimental Realizations 

Ignoring the nonlinear term, one can perform the sums in Q5.187D and fl5.18SQ numerically to 
find the magnetic field and lattice structure. The correction due to the nonlinear term can then 
be found iteratively. To do so, we start by writing ( 5.184| ) in Fourier transformed form (with 
Gaussian source) as 

R n e-^ + eQ^- 1 ] 

where Q[B&] is the Fourier transform of Q[B] and the superscript n represents the stage of 
iteration. The initial B9 is found by neglecting the nonlinear term. At each following step we 
calculate Q[Bk] using the known value for Bk, substitute it into (5.189) to find the new B^. We 



continue this iteration until we get the desired accuracy. Having B, we then can calculate the 
Gibbs free energy Ql = Fl — HB / 4ir using (|5.183| ). We calculate the Free energy this way for 
different configurations of the vortex lattice. We find that a centered rectangular flux lattice, 
with principal axes aligned with the ionic crystal lattice minimizes the free energy. An angle 
P between unit vectors (Fig. 5.10| ) characterizes the dependence of the lattice on e and the 



magnetic field. An example of such a centered rectangular lattice is shown in Fig. 5.1C 



In agreement with earlier results within GL [121] and Eilenberger |147j] formalisms, indi- 



vidual vortices are elongated along the crystalline axes. Figure 5.11| (a) shows the dependence 



of Gibbs free energy on j3 for various values of e at fixed applied field H = 400-Bo — 6.8T. For 
e = minimum occurs for /3min = 60°, corresponding to a hexagonal lattice. As e increases, 
/?min continuously increases and for sufficiently large e, the flux lattice becomes tetragonal with 
/?min = 90°. For /?min 7^ 90°, there are always two solutions, related by a 90° rotation, in which 
the long axis of the centered rectangle is aligned with either the x or y axis. The degeneracy is 
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Figure 5.10: Distribution of magnetic field in a vortex lattice for e = 0.3 and H = 6.8T, 
leading to an angle of /3 ~ 74°. We use Ao = 1400A and k = Ao/£ = 68. 

much larger for e = 0, when the flux lattice may have an arbitrary orientation relative to the 
ionic crystal lattice. 

The dependence of /?min on the applied field for various values of e is displayed in Fig. 
5.1 l| (b) . Clearly the anisotropic term becomes more important at larger fields. Our perturbative 



elimination of v in favor of B breaks down when e and H are sufficiently large that /?min differs 
significantly from 60°. Furthermore, we might expect higher order corrections to fl5.177| ) to be 
important in this regime. By fitting Fig. 5.11| (b) to experimental data on tetragonal materials 



such as Tl2Ba2CuOg +( i (once such data become available) one can directly assess the magnitude 
of e, the only unknown parameter in the model. 

The analysis presented here, can be easily extended to take into account effective mass (i.e. 
penetration depth) anisotropy. In a simple one-component GL model, the derivative term is 
generalized to: 

f= 7 4 |IM 2 . (5.190) 

i=x,y,z 

We restrict our attention to fields along the z-axis. The anisotropy then can be removed by a 
rescaling of the x-coordinate and a corresponding rescaling of the magnetic field. The coherence 
length and penetration depth anisotropies are considered the same: £ y /£, x = A x /A y . Making 
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Figure 5.11: a) Gibbs free energy as a function of (3 for the same parameters as Fig. [5.10| 
and various values of e. Arrows indicate positions /3min of the minima and Go = —H 2 /87r. b) 
Equilibrium angle /3min as a function of H for several values of e. 

a simplifying assumption that the higher derivative and mixed derivative terms in J- are also 
simply modified by a rescaling by a common factor, it then follows that the flux lattice shape is 
obtained by stretching along the x-axis by the factor \ x /\ y . We now obtain two possible vortex 
lattices, both of centered rectangular symmetry, aligned with the ionic lattice, with different 
angles, j3. (Relaxing our simplifying assumption may split the degeneracy between these two 
lattices.) On the other hand, when e = 0, we may rotate the hexagonal lattice by an arbitrary 
angle before stretching. This gives an infinite set of oblique lattices with arbitrary orientation. 

To compare theory with YBCO we should take into account twin boundaries; which are the 
boundaries between different orientations of crystal lattice. They may also tend to align the 
vortex lattice by pinning vortices to the twin boundaries, at ±45° to the x-axis. This effect 
competes with alignment to the ionic lattice which we have been discussing. Only in the special 
case of a square vortex lattice does a line of vortices occur at ±45°. If this is not the case, 
and if pinning by twin boundaries is significant, then we should expect that the vortex lattice 
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will align with the ionic lattice far from twin boundaries but will be deformed in the vicinity 
of a twin boundary in an effort to align itself with the twin boundary. On the other hand, 
for e = 0, the vortex lattice would remain aligned with the twin boundaries everywhere except 
within vortex lattice domain boundaries which necessarily exist roughly midway between the 
twin boundaries. 



Neutron scattering experiments on YBCO [117] suggest that the vortex lattice is well-aligned 
with the twin boundaries and is close to being centered rectangular (the ratio of lattice constants 
is about 1.04) with (5 ~ 73°, with weak dependence on H. This corresponds to a rotation away 
from alignment with the ionic lattice by 9°. Four different orientational domains, related by 
reflection in the (1, 1,0) axis and 90° rotation were reported. These results can be rather well 
fitted [124] by the basic London model (e = 0) with mass anisotropy. For X x /X y = 1.5, a value 
roughly consistent with infrared experiments, this lattice has about the right shape. Taking 
into account the two crystallographic domains (related by interchanging and X y ) there are 
altogether four vortex lattice domains, as seen experimentally. The experimental fact that the 
vortex lattice appears to be well aligned with the twin boundaries suggests that the tendency 
to align with the ionic lattice is small. No evidence for a bending of the vortex lattice (by 9°) 
into alignment with the ionic lattice far from the twin boundaries has so far been found. 



STM imaging of the YBCO vortex lattice | 118 | also suggests that the (highly disordered) 
lattice has approximately centered rectangular symmetry with (5 ~ 77°. However, no evidence 
for the 9° tilt into alignment with the twin boundaries was reported. Considering the observed 
anisotropy of the vortex cores it has been concluded that the mass anisotropy alone cannot 
account for the measured 77° angle of the vortex lattice. It has been suggested that a mechanism 
related to the internal symmetry of the order parameter (such as the one discussed in the present 
thesis) needs to be invoked in order to reconcile these observations. 



Low field Bitter decoration data on YBCO [148] show vortex lattice geometry with a very 
small distortion from hexagonal, consistent with a much smaller anisotropy \ X /X y = 1.11 — 1.15. 
One may be tempted to attribute this apparent field dependence of (3 to the effects discussed 
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above in connection with Fig. 5.11| (b). An alternative explanation is a poor quality of samples 
used in the Bitter decoration experiments that may have resulted in partial washing out of the 



a-b plane anisotropy [14£]. 
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Figure 5.12: Magnetic field distribution function; Solid line: using the same parameters as in 



Fig. 5.10. Dashed line: the same parameters except for e = 0. 



/iSR experiments measure the field distribution n(B) = (l/f2) / 5[B — B(fj\dxdy. Our 
calculation of this quantity is shown in Fig. [5.12 . For comparison we show our results with 
(solid line) and without (dashed line) the extra nonlocal term. For (3 ^ 60°, B(r) has two 
inequivalent saddle points leading to two peaks in n(B). n{B) is unaffected by effective mass 
anisotropy, as can be shown by the rescaling transformation, mentioned above. Existing fiSR 



experiments show only a single peak [15C] but this might be due to the large broadening as a 
result of finite muon lifetime. 

The weak field dependence of /3, the alignment with twin boundaries in the neutron scat- 
tering experiments, and the single peak in n(B) suggest that e is small in YBCO and that the 
normal London model, together with twin boundary pinning, provides a good fit to the data. 
STM and Bitter decoration data on the other hand, seem to favor finite e and weak pining to 
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twin boundaries. Further experimental work, preferably on untwinned YBCO or other tetrag- 
onal superconductors, will probably be necessary to clarify the importance of square lattice 
anisotropy in high-T c superconductors. 

Most of the experiments mentioned above are done at temperatures far below T c , i.e. beyond 
the validity range of GL theory. Thus our derivation of the generalized London model from GL 
theory is not valid for most experimental situations. However, as we emphasized, the corrections 
we introduce to the London free energy are the most general analytic corrections possible by 
symmetry to the order of our concern. It is therefore reasonable that even at lower temperatures 
all the lowest order analytic corrections come with these forms. We will see in the next chapter 
that at very low temperatures, the lowest order corrections are actually non-analytic. This is 
because of the existence of the gap nodes; a unique feature of d-wave superconductors. Even 
those effects turn into the analytic forms discussed here at higher temperatures. 



Chapter 6 



Microscopic Model 



6.1 Gorkov Theory and Quasi-Classical Approximation 

Ginsburg Landau theory, discussed in the last chapter, is a theory that is valid only near the 
transition temperature T c . Most of the experiments however, are done at temperatures far below 
T c . Also, one of the most essential characteristics of d-wave superconductors, i.e. the existence 
of nodes in the superconducting gap, does not show a direct effect on the GL equations. Thus 
a more sophisticated microscopic approach is necessary to describe the vortex lattice properties 
of high T c superconductors. The microscopic theory we present here is based on the theory 



developed by Gorkov [151] in 1959. In his paper, Gorkov showed how GL theory results from 
pairing theory at temperatures close to T c . Here, we don't give the detailed derivation of 
Gorkov's equations from the first principles. Our goal instead, is to derive relations between 
supercurrent j and superfluid velocity v s (or vector potential A) starting from these equations. 
This relation is then used to obtain a generalized London equation. For a detailed derivation 
of the Gorkov equations, interested reader should refer to ref. [136, 152). 
Let us start from an interaction Hamiltonian of the form 

H int = J d 3 rd 3 r'^\{^)^{f^)^\{r)'ip l {r)V{r-r') (6.191) 

The Gorkov equations for this Hamiltonian are 

+ J A(fy')^{7'TyT')d 3 r" = 5 {3) (r - t*)8(t - r') (6.192) 

{^ + \^ + f^)] 2 + M}^(rr, fV) 

A*(r,r")g(r"T,r<T')d 3 r" = (6.193) 
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where Q and T are regular and anomalous parts of the Nambu-Gorkov Green function respec- 
tively, 

S(fT,r-V)<W = - <T r Vv(rr)Vi,(rV) > (6.194) 
T\rT,rT) = - <T T ^{(rr)^!(rV) > (6.195) 

Here, and also throughout the rest of this chapter, we choose h = m = 1. The order parameter 
A is related to the interaction potential V by 

A(f - = V{r- f) < V>t WltrO > • (6.196) 

Equations (|6. 192f) and ( |6.193| ) have to be solved for Q and J- self-consistently with A. The 
self-consistency equation follows from ( |6.195| ) and ( |6.196 ) to be 



= V{r-r')^{r0 + ,r f 0) (6.197) 

Solving these equations self-consistently, in a general case at the presence of Magnetic field 
is tedious if not impossible. However, substantial simplification can be achieved using local 
approximation or gradient expansion. The method we use here (with some minor differences), 
is a generalization of the gradient expansion method discussed in Ref. |152(| . The generalization 
is necessary to study d-wave superconductivity. 

6.1.1 Local Approximation 



Gorkov equations ( 6,192| ) and ( |6.193| ) take much simpler forms in Fourier space. If the system 



is translationally invariant in time, Q and J- are only functions of (r — r'). We can therefore 
take their Fourier transform with respect to (r — r') and define 



GUr,r) = J G{rT,70)e- WnT dT 

?l(?J) = J P(rT^0)e- iuJnT dT (6.198) 
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where io n = nT(2n — 1) are fermionic Matsubara frequencies. We also change the coordinates 
to R = r' and p = r — r' and take the Fourier transform with respect to p to get 

g u (R, k) = j g u (R + p, R)e~ il Pd 3 p (6.199) 
PjR, k) = J Tt{R + p, R)e~ il Pd 3 p (6.200) 

A(R, k) = j A(R + ?R- ^e'^dPp (6.201) 
J 2 2 



( 6.192 ) therefore becomes 

1 r.r ie .8 



{iuj n + hik--A(R + i^)} 2 + p}gUR,k)+ f d 3 p f d 3 r"e- il PA(R + p,r")^' 1 R) = l 
2 c dk J J 

(6.202) 

The integral on the left hand side of ( |6.202 ) can be written as 



J d 3 Pl J d 3 p 2 e- il ^+^A(R + p 1 + p 2 ,R + p 2 )Fi(R + p 2 , R) 

= lim / tfVi / d 3 p 2 e- Ck '^e~^ A((R + p 2 + ^) + (R + p 2 + ^) - ^JR + p 2 , R) 

k'-*kJ J I I I I 

= Yim^A{R + i^ + ^k')Tl{R,k) (6.203) 
k'^k ok 1 ok! 

If A(R) and A(R, k) are slowly varying functions of R, then we can neglect all and which 
are at most of the order of £. We can also follow the same procedure for ( |6.193| ). The resulting 
Gorkov equations in this limit will be 

{iu n -hk--A(R)] 2 + p}gUR,k) + A(R,k)^(R,k) = 1 (6.204) 
z c 

{-iuj n -hk + -A(R)] 2 + p}^(R,k)- A*(R,k)g u (R,k) = (6.205) 

These are now algebraic equations instead of integral equations. This approximation is usually 
referred to as quasi-classical approximation. The reason is that the equations involve position R 
and momentum k at the same time, similar to classical equations. The self-consistency equation 
also takes the form 

A* (R, k) = i Yl I d 3 k'V(k - %)J* (R, k') (6.206) 

P n. J 
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where V[k) is the Fourier transform of V(r). The solutions to ( |6.204j ) and ( |6.205| ) are 



-iu n -\[k + %A{R)f + n 



{iuj n - - f A{R)Y + n}{-iuj n -\[k+ ZA(R)} 2 + n} + A 2 (k) 

Tl{R,k) 



{ iuJn _ i[ k _ |A(i?)] 2 + - + fl(i?)] 2 + m} + A 2 (fc) 

In the absence of magnetic field these solutions will be simplified to 

—ioj n — e? 

k 

= (6-207) 

k 

with 



H = X -k 2 -ii and £g = Je| + A 2 (fc) (6.208) 
The self-consistency equation ( |6.206| ) also will be 

A*(R,k) = [ d 3 k'V(k-k')A*(k')W^TW (6 - 209) 
Summing over Matsubara frequencies, the self-consistency equation becomes 

A(k) = -J d 3 k>V(k - k')^- tanh(^) (6.210) 

which is the well known BCS gap equation. 

In the presence of magnetic field, simplifications can be achieved if we introduce new Green's 
functions Q and T by 

g(f r ,fV) = Q(fT^T l )e- C(x ^ ) - x ^ )) 

^(fr,fV) = ^(fy,rV)e i(x(f)+ * (<,,)) (6.211) 

and also a new pairing potential 

A(r, /) = A(r, F) e ^ l (x(0+x(^)) ( 6 2 i 2 ) 
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where xi?) ls the phase of the operator ip(r) at some fixed gauge. The Gorkov equations fl6.20-| ) 
and ( 6.205| ) for Q and J- become 



{^ n -i[fc + t; s (i?)] 2 + ^(ii,fc) + A(i?,fc)^t(i|^) = i 
{-iu; n -hk-v s (R)} 2 + (i}ft(R,k)-A*(R,k)GUR,k) = 



(6.213) 
(6.214) 



where v s is the superfluid velocity defined by 



6 -* 

v s (r) = Vx(r) Air) 

c 



(6.215) 



This is the same definition as in ( 5.148 ) if we notice that x = 0/2- Shifting the chemical 
potential // by Vg/2, we have 



1 



(k ± v s ) 2 - \i = eg± h- v s 



Gorkov equations will then be simplified to 
/ 



V 



and their solutions will be 



iLO n -e^-k- v s 
A* 



iuj n + er - k ■ v & 



G 



(6.216) 



(6.217) 



Gu>{R,k) 



iuj n - k- v s (R) + 
[iuo n -k-v s {R)f -El{R) 

-A*(R,k) 
[iu n -k-v s (R)f -El(R) 



(6.218) 



Notice that the only difference between these and (|6.207|) is a Doppler shift (by k ■ v s ) in the 
Matsubara frequencies. We again sum over the Matsubara frequencies in the gap equation to 
get 

A(R,k) = - [ d 3 k'V(k - k') A{K fj {f[k' • v s (R) - E p (R)} - f[k' • v s (R) + E p (R)}} (6.219) 



where 



f(E) 



1 



(6.220) 



is the Fermi distribution function and (3 is the inverse temperature 1/T. Note that as v s — > 0, 
this reduces to ( |6.210| ). 
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6.1.2 Calculation of Current in a Superconductor 

The current density, in terms of the Green functions, can be written as 

o 

— > i e x p — * x 

j(r) = - lirn £(V^ - V P ) < ^(^ > --A(f) ]T < ^(f)^(f) > 



c 

2 

lim ^(V F - V f )g(rD, f'0+) - —A(r) ]T <5(rt), rD+) 



2 r'—>r' 

u u 

= ie lim(V r -/ - V r -)a(H),r + ) + 2ev s (f)g(rO,rV + ) (6.221) 

r— >r* 

The current in this form is manifestly gauge invariant. Keeping in mind that the superfluid 
density is n s (r) = < r4-(j)Vv(^) >= 2Q(fO,rO + ), the current can be written as 

j{r) = en s v s {r) + j qp (r) (6.222) 

where 

j (f) = ie lim (Vf» - V r -)£(rt), ^0+) (6.223) 

r— >r' 

j qp turns out to be the current due to quasiparticle excitations in the superconductor. Fourier 
transforming, using the same techniques discussed in the last subsection, we get 



-> f d?k 1 -> — 

Jqp{R) =ie J J^ ( - V A- 2ik hT,^(R,k) (6.224) 



(2vr) 3V R 'fi 

Again under quasi-classical approximation, one can neglect and write 
2e^ f d 3 k 



j qP (R) = 7 E/^iM*) 

/rfik -> ^ ^ 

--3fc{(l + f)f[E- k + fc • v s {R)\ - (1 - ^)/[^ - k ■ v s (R)}} (6.225) 



In the last step, we have performed a Matsubara sum. Combining the two terms the current 
can be written as 

- f d 3 k - 

j qp (R) = -2e J j^kf[E S - k ■ v s (R)] (6.226) 

or 

j qp (R) = -AeN F (v F J™ dU (yWlAP - v F • t/ a ^ (6.227) 
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where Np is the density of states at the Fermi surface and (...) FS represents averaging over 
the Fermi surface. Note that A, in general, can depend on the point on the Fermi surface. 
Equation ( |6.227| ) has a natural interpretation. The argument of the Fermi distribution function 
/ is the Doppler shifted quasiparticle energy at the particular point on the Fermi surface. Thus 



(6.227) is actually a Fermi surface average of the current due to the quasiparticles excited via 
thermal excitation plus the Doppler shift. This can be viewed as a normal fluid in the opposite 
direction of the superfluid that tends to reduce the supercurrent density. At zero temperature, 
the Fermi function is a step function which is zero unless its argument is negative. Therefore 
d6.227f) is zero unless 

v F -v s > |A| (6.228) 

Eq. (|6.228; ) actually defines a critical value v c = |A| m i n /u^ for the superfluid velocity v s , 
beyond which the superconducting gap vanishes at some regions of the Fermi surface. Increasing 
v s beyond v c increases the number of quasiparticle excitations and eventually destroys the 
superconductivity. Therefore, there is an interval of v s over which we can have superconductivity 
with a gap that vanishes in some region of the Fermi surface. This is usually known as gapless 
superconductivity. For an s-wave superconductor, Q6.228 ) gives a finite critical value for the 
superfluid velocity. For a d-wave superconductor (or any superconductor with nodes in the 
superconducting gap) on the other hand, the gap vanishes at the nodes even at zero v s . In 
other words, quasiparticles exist at node lines even for infinitesimal v s and zero temperature. 
This phenomenon is highly anisotropic because the nodes exist only at four symmetric points 
on the Fermi surface. The effect of these nodal quasiparticles in the Meissner state was studied 



by Yip and Sauls [104] and others ||153|j. Here, we focus on the effect in the vortex state; which 



has been studied in [132, 146]. We again derive a generalized London equation and solve it 
numerically for a vortex lattice to calculate physical quantities. 
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6.2 Generalized London Equation 

As we mentioned earlier, the fourfold anisotropic nature of the gap nodes is expected to influence 
the magnetic properties of a d-wave superconductor. This influence will be via two distinct 
effects: a nonlinear effect due to quasiparticle generation because of Doppler shift in the quasi- 
particle spectrum, and nonlocal effect due to divergence of the coherence length along the nodal 
directions. In this section we study these effects in two separate subsections. In what follows 
we assume that these two effects can be studied independently and their effects are additive. 
In other words, we assume that any inter-connection between these two effects produces higher 
order corrections to our generalized London equation. In general the excitation of quasiparticles 
at the gap nodes can affect the coherence length and therefore the nonlocality of the system. 
Similarly nonlocalities can change the form of Eq. fl6.227| ). This effect which is neglected in 



our calculations has been studied extensively in the context of superfluid He [154]. 



6.2.1 Nonlinear Corrections 

Let us first neglect any nonlocal effect and focus on the quasiparticles generated at the gap 
nodes. As we mentioned at the end of the last section, excitation of quasiparticles at the gap 
nodes produces a current density flowing in the direction opposite to the superfluid velocity - 
sometimes called back-flow. The total current is given by fl6~2"22l ). We again define the superfluid 
velocity by 

v s = -(V(/}-— A). (6.229) 
2 c 

with <p being the phase of the order parameter. In this definition, v s is in the direction of 
the superfluid, i.e. in the opposite direction of motion of the Cooper pairs (which is used 
in |104| , |146j ] as the definition for the direction of v s ). The contribution of the quasiparticles 
generated at the nodes to the total current is given by (|6.227 ) which we rewrite it as 



j gp = -AeN F f dsv F (s) I" dU{Je + Hs) 2 ~ M*) ' 30 ( 6 - 230 ) 
Jfs Jo v 
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Figure 6.13: Circular Fermi surface with a d x 2_ y 2 gap. Quasiparticles will be excited at nodes 
marked by v F \ and vp2 opposite to v s . 

where s parameterizes a point on the Fermi surface. A(s) is the superconducting gap which in 
general can have s-wave, d-wave, or other symmetries. At zero temperature ( |6.230 ) leads to 



j qp = -4eN F J dsv F (s)9{v F ■ v s - | A|) J (v F ■ v s ) 2 - A 2 (6.231) 

At higher temperatures however, this gives the first term in the Sommerfeld expansion; 
which is a good approximation as long as T < T* = T c (H/Hq), where Hq is of order of the 



thermodynamic critical field, H c 104]. The presence of the ^-function in ( 6.231 ), results in 



excitations only at the nodes which are in the direction of v s . Fig. 6.12 illustrates a circular 
Fermi surface with a d-wave gap. The quasiparticles are excited at the nodes marked by v F \ 
and v F 2 in the direction of v s . For a small enough v s , the excitations stay very close to the gap 
nodes. Therefore, one can linearize the gap function near the nodes writing A(#) ~ 7Ao#, with 
Ao the maximum gap and 7 defined by 

7 = ^[^A(^)] node (6.232) 

The most commonly used form for a d-wave gap is A(0) = Aocos(2#). In that case, ( |6.232| ) 
leads to 7 = 2. The component of j qp along the x'-direction which is diagonal to the x and y 
(a and b) directions (as illustrated in Fig. 6.13| ) is then 



f 6c d6 I 

jqpx' = -4eN F v F J — \j {v F v s cosa) 2 - \^Aq9\ 2 
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-en s — cos a cos a 



-en. 



^sx' I ^sx' 



(6.233) 



a is the angle between v s and x'-axis, n s = Npv 2 F is the superfluid density, vq = jAq/vf is 



some characteristic velocity and 9 C is a cut off imposed by the ^-function in ( |6,231 ). Similarly 
the ^/'-component is 



Jqpy' 



-en. 



and the total current thereby becomes 



en s [v s - {v sx >\v sx >\x + v sy '\v sy >\y')/vo\. 



(6.234) 



(6.235) 



The nonanalytic nature of the effect is evident from this equation. It is now possible to write 
a free energy in such a way that ( gggp can be obtained by minimization with respect to 
A. Keeping in mind that j = (c/4ir)V x B = (c/47r)V x V x A and d/dv s = —~d/dA, the 
corresponding free energy density can be written as 



^ s 2 -^(k*'l 3 + K' 



+ 



8vr 



(6.236) 



In general, it is possible to solve ( |6.235 ) for v s in terms of j = (c/47r)V x B, substitute it into 



(6.236) and write down a London free energy only in terms of B and its derivatives. However, 
instead of solving ( 6.235j) exactly, we find v s perturbatively assuming that the nonlinear part 
is much smaller than the linear part. This way we get a polynomial correction to the London 
equation which is convenient for numerical purposes. To first order in perturbation theory we 
have 



47reno 



V x B + 



(V x B) X >\V x B\ x ,x' + (V x B) y t\V x B\ y >y'] 



(6.237) 



4:7ren s vo 

Substituting this into ( |6.236| ) and keeping the lowest order terms, the London free energy density 
becomes 



h = ± [B 2 + Ag(V x Bf + (|^(I(V x B) x ,f + |(V x B) y ,f) 



(6.238) 
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where Ao = \/c 2 /47re 2 n s is the zeroth order penetration depth, £o = vf/ttAq is the coherence 
length, Bq = (/>o/2-7rAo is a characteristic field of the order of H c \ and <po = irc/e is the flux 
quantum. For magnetic fields in the z-direction, ( |6.238 ) becomes 



h = ±- \B 2 + X 2 (VBf + A^(\d x ,Bf + \d y ,Bf)} =f L + U (6.239) 

07T L 37 Bq J 

with /£ and / n ; representing the ordinary London free energy density and the leading nonlinear 
correction to the free energy density respectively. The corresponding London equation is given 
by 

- \ 2 X7 2 B + B- (—)^(d 2 x ,B\d x ,B\ + d 2 v ,B\dy,B\) = (6.240) 

J B 

A similar London equation is also derived by Zutic and Vails who investigated the effect in the 



Meissner state |153|]. 

In order to find the magnetic field distribution in a vortex lattice, one has to insert a source 
term Y>jp(f — fj) on the right hand side of ( |6.240| ) with fj being the position of the vortices in 
the lattice. The function p(r) takes into account the vanishing of the order parameter at the 
center of the vortex cores. Numerically, it is more convenient to work in Fourier space rather 
than real space. Fourier transforming ( |6.240| ) with a proper source term on the right hand side 
yields 

B % + A^fc 2 ^ - G nl (k, B % ) = BF(k) (6.241) 

where A; is a reciprocal lattice wave vector, G n \ is the Fourier transform of the nonlinear term and 
B is the average magnetic field. The cut off function F(k) comes from the Fourier transformation 
of the source term and removes the divergences by cutting off the momentum sums. We are 
going to use Gaussian cutoff ( |5.168 ) in our numerical calculation in the next section. 



6.2.2 Nonlocal Corrections 

In our derivation of the corrections to the London equation in the previous subsection, we 
derived a local relation between current density j and superfluid velocity v s (or vector potential 
A). This was a result of the local approximation we used in our derivation. In fact, this relation 
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is always nonlocal over the length scale of the coherence length £oj which is of the order of the 



finite spatial extent of the Cooper pair 130|. Magnetic field in a superconductor varies in the 
length scale given by London penetration depth Ao and therefore nonlocal corrections to physical 
quantities, such as the effective penetration depth, will be of order k~ 2 , where k = Ao/£o is the 
GL ratio. For strongly type II materials (k 3> 1) such corrections are negligible. Since cuprate 
superconductors fall well within this class (k is of 50 for most) local electrodynamics is usually 
used. However, a closer examination suggests that this might not be justified in all situations, 
if, as it is widely believed, these materials exhibit nodes in the gap. In such a case in place of 
the usual coherence length £o = vf/kAq one is forced to define an angle dependent quantity, 
£o(j5) = vf/ttAp, which diverges along the nodes. Clearly, in the vicinity of nodes the condition 
\o/£,o(fi) 3* 1 is no longer satisfied and, in fact, the extreme nonlocal limit is achieved. Nonlocal 
corrections therefore cannot be dismissed in unconventional superconductors. The effect of 
these corrections on the temperature dependence of penetration depth at low T was studies by 
Kosztin and Leggett [111]. Here, we discuss about the effect on the vortex lattice properties. 



What we follow here is based on a work published in Ref. [132]. From the above arguments, it is 
clear that corrections due to the nonlocal effects at the gap nodes will be highly anisotropic and 
will in general break the rotational symmetry of the flow field around the vortex, contributing 
an anisotropic component to the inter-vortex interaction in the mixed state. 

Let us now find these nonlocal corrections to the London model neglecting all the nonlinear 
effects discussed in the previous subsection. In the next subsection, we will combine these two 
effects into a single London equation to be used in numerical evaluations. Nonlocal relation 



between j and A is conveniently written in the Fourier space as [ |13C| 

j k = -(c/4n)Q(k)A k . (6.242) 

Here Q(k) is the electromagnetic response tensor which can be computed, within the weak 
coupling theory, by generalizing the standard linear response treatment of Gorkov equations 
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[155] to an anisotropic gap. We find 



Q*® = ^ E / ; *** Fi * Fi ) , (6-243) 

3 A o ^o\ v A4Ta|k + a| +7 |)/' 



where 'y^ = vf ■ k/2, Ao = <J c 2 / ' ^Ke 2 v 2 F Np is the London penetration depth as before, u> n are 
the Matsubara frequencies and the angular brackets again mean the Fermi surface averaging. 



Derivation of ( 6.243 ) involves a calculation of the perturbative correction to the Green's func- 
tions to linear order in A using Gorkov's equations. The result is then used to calculate the 
supercurrent j |155| ]. Eq. ( |6.243| ) is valid for arbitrary Fermi surface and gap function. For 
isotropic gap one recovers an expression derived by Kogan et al. from the Eilenberger theory 
|156f| . One may simplify solving the London equation again by writing it in terms of magnetic 



field only. Eliminating j from (6.242) using the Ampere's law j = (c/47r)V x B, one obtains 



B k - k x [Q~ l (k)(k x B k )] = 0. (6.244) 

To study vortex lattice using this formalism, it is necessary to insert a cutoff function F{k) on 
the right-hand side of the Eq. ( |6.244[ ). In our numerical calculations in the next section we 
again use Gaussian cutoff function ( |5.168| ). For many purposes it is also convenient to write 
down the corresponding London free energy, such that bThj^Bk = gives the above London 
equation: 

T L = Y\B\ + (kx B k )Q-\k)(k x S fc )]/87T. (6.245) 

k 

At long wavelengths Q(k) can be evaluated by expanding expression ( |6.243| ) in powers of 
7S. The zeroth order term 

k 

Q " =s " x ^Sxm+asjW' <6 ' 246) 

clearly recovers the ordinary London free energy. Here, A = A(T) is the temperature dependent 
penetration depth which, at low temperatures, has the well known T-linear behavior [157, 33]; 



for a d x 2_ y 2 superconductor with = Ad(p 2 — p 2 ) and a cylindrical Fermi surface (from now 
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on we shall focus on this simple case) A 2 ~ A 2 (1 — 2 In 2T/A^). The leading nonlocal term is 
quadratic in k: 



The expression Qij = 5ij\ 2 + is easily inverted to leading order in k to get « 
A 2 [<% — X 2 q\^}. Substituting this into ( |6.245| ) and specializing to fields along the z-direction 



A 2 1 r 27r A 2 w,, 



we have 

T L = B \[ l + + A 2 £ 2 (ciA: 4 + c 2 k 2 x k 2 )]/8ir. (6.248) 

k 

Here £ = vp/nAd and A^ is assumed to be a temperature dependent solution to the appropriate 
gap equation. Dimensionless coefficients c\ and C2 are given by 

\2 1 r 2w A 2 w, 

n>0 

where w\ = v 2 Fx v 2 Fy , W2 = (v Fx — v Fy ) 2 — ^v 2 Fx v Fy and the Fermi surface has been explicitly 
parameterized by the angle 8 between p and x-axis; vf = (cos 8, sin 8) and A^ = A^ cos 28. 
Coefficients c\ and C2 depend on temperature through a dimensionless parameter t = T/Aj. 
From ( UH) one can deduce their leading behavior in the two limiting cases: for iClwe find 

7T 2 A_2_1 

and for t 3> 1 (i.e., near T c ) 



c i = ^v27 « c 2 = "4ci, (6.250) 

O An t 



A 2 1 

Cl = a ^2^4 ' c 2 = 8ci, (6.251) 

where a = C(5)(l — 2 _5 )/87r 2 = 0.01272. In the above A also depends on t, but this will be 
unimportant for the following qualitative discussion. 

The free energy ( 6.248| ) formally coincides with ( 5.188j ), deduced previously in the last chap- 
ter with phenomenological considerations. The present model however, allows us to evaluate 
the coefficient of the symmetry breaking term C2 (or e in ( [5.188] )). At high temperatures, C2 is 
found to be positive; in agreement with our previous result. As we showed in the last chapter, 
such a term leads to a centered rectangular vortex lattice structure with principal axes oriented 
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Figure 6.14: Two high symmetry orientations of centered rectangular unit cell. 



along x or y axes of the ionic lattice (see Fig. |6.14| a). The magnitude of distortion from a per- 



fect triangular lattice is controlled by the magnitude of C2 and grows with increasing magnetic 
field. Eq. ( |6.251| ) shows that at fixed field this distortion will initially grow with decreasing 
temperature. At low temperatures ( |6.250 ) predicts C2 < 0. This will lead to the same centered 



rectangular lattice but rotated by 45° (Fig. |6.14 b). Numerical evaluation of ( |6.249| ) shows that 



C2 passes through zero at if ~ 0.19. At this temperature the free energy (6.248) is isotropic 
and the lattice will be triangular at all fields. The change of sign of C2 reflects the competition 
between the two terms of different symmetry in W2. At t <C 1 only the region around the node 
is important where (v Fx — Vp y ) 2 vanishes, while at t 3> 1 the average involves the entire Fermi 
surface to which both terms in W2 contribute. This change of sign is a unique consequence of 
nodes in the gap function and would not occur in conventional superconductors. 

Another consequence of nodes is the fact that, as can be seen from ( |6.250D , both c\ and C2 
diverge as 1/t for t — * 0. This divergence signals that the response tensor Q(k) is a non-analytic 
function of k at T = 0, and the expansion in powers of 7S breaks down. Formally this is caused 
by the fact that at T = 0, at the nodal point, the expression ( |6.243| ) for Qij(k) contains a term 
proportional to l/7§- At T = the frequency sum in ( |6.243j ) becomes an integral which can be 
evaluated exactly with the result: 

~ /7\ 1 / ~ „ 2arcsinhw\ 

Qij(k) = ^2 ( v Fi v Fj / ) , (6.252) 

A o \ yVl + y 2 / 
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where y = j^/Ap. For small k the dominant contribution to the angular average comes from 
the close vicinity of nodes and can be evaluated by linearizing around the nodes. One 
finds that the leading nonlocal contribution is linear in k rather than quadratic. For Qij = 
<%^0 2 + Qif> we nave Qxx = Qvy = - M(^mCo) and Qxy = Q y l x = -n(k m ^ )sgn(k x k y ), where 
k^i = max (| |, \k y \) and k m = mm(\k x \, \k y \). Prefactor \i = 7r 2 /8\/2 = 0.8723 is exact in the 
sense that all corrections to Qij are 0(k 2 ). The resulting free energy at T = is 

^ = E + A o fc2 + Mo^M^M - kl)]/8ir. (6.253) 

The nonlocal term is clearly non- analytic in k. Its functional form is universal in the sense that 
it is independent of the Fermi surface structure (as long as it has the same symmetry as the 
order parameter) and the prefactor /i only depends on the angular slope of the gap function 
and Fermi velocity at the node. 

The present calculation can be easily generalized to treat the effects of Fermi surface 
anisotropy. As mentioned above tetragonal anisotropy will not modify the T — ► universal be- 
havior but may lead to quantitative changes at higher temperatures. Orthorhombic anisotropy, 
on the other hand, will modify even the T — * limit. We expect that it will, to leading order, 
merely rescale the coordinate axes, leading to the same structures as described above stretched 
by the axes, leading to the same structures as described above stretched by the appropriate 
factor [ p-40| . It may further remove the degeneracy between two equivalent lattices related by 
90° rotation. Also neglected in our calculation is the effect of electronic disorder, which will 
remove the non-analyticity of Q(k) at longest wavelengths, just as small finite temperature 
would. Since the lattice structure is most sensitive to Q(k) at finite k ~ l/l (I is the vortex 
spacing), we expect our predictions to be robust with respect to weak disorder. 

6.2.3 Combining Nonlinear and Nonlocal Corrections 

Having established theories for nonlinear and nonlocal effects independently, we now would like 
to combine them in one London equation to be used in numerical calculations. Nonlocal effect 



Chapter 6. Microscopic Model 



95 



can be added to the generalized London equation (6.241), simply by replacing the second term 
with its nonlocal counterpart 



Br + dji^kikjBr - G nl (k, Br) = BF(k). 



(6.254) 



Sums over i and j are implicit here. £ij(k) is related to the electromagnetic response tensor 
Qik) by 

C i:j (k) = Qij{k)/ det Q(k). (6.255) 

The electromagnetic response tensor Qik) is given by ( 6.243 ) for T / 0, and by ( 6.252| ) at 
T = 0. 

In the local case, we have £ij(k) = which leads back to ( |6.241| ). Taking G n \ to the 



right hand side of (|6.254[) , Br can be obtained by 



B T 



BF(k) + G nl (k,Br 



(6.256) 



1 -\- Hij(k)k{kj 

We will use this equation, in the next section, to find B^ iteratively for a specific lattice geometry. 
Having B^, the free energy can be easily calculated using 



F = ^nl + J} 1 + Ciji^kik^B 



(6.257) 



where J- n \ is the free energy due to the nonlinear term in Eq. ( S.239j ) 



3 7 ' B 



d-r (\d xl Bf + \d y ,B\ 



(6.258) 



6.3 Numerical Calculations 



Unlike Ginsburg-Landau, the London free energy cannot completely determine the vortex lattice 
by a simple minimization. Instead, one has to impose a set of source terms located at the 
position of the vortices on an assumed lattice. The functional form of the source terms does 
not come from the London theory and requires more fundamental treatments; as we discussed 
in Sec. (5. 2. 2). Information about the functional form of the source term is included in ( |6.256| ) 
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via F(k), and information about the positions of the vortices is included in the reciprocal lattice 
vectors k. The free energy must then be minimized with respect to the positions of the vortices 
in order to find the equilibrium lattice configuration. In general, a 2D lattice can be determined 
by four parameters. However, it turns out that a centered rectangular lattice is energetically 
more favorable than an oblique lattice. On the other hand, the vortex lattice spacing is fixed 
by the average magnetic field B (B ~ H away from H c \). Thus we are left with two variational 
parameters, i.e. the lattice orientation with respect to a and b directions and the apex angle 
(3 - the angle between the two basic vectors of the lattice. We therefore find the vortex lattice 
geometry by minimizing T in with respect to the apex angle (3 for different orientations of the 
lattice. 

In our numerical calculation, we first neglect the nonlinear effects and consider only the 
nonlocal corrections. At the end, we will add the effect of the nonlinear corrections to the results. 
We shall see that the nonlinear corrections are subdominant to the nonlocal ones. Therefore 
the results we obtain, ignoring the nonlinear effects, are actually very close to the reality. In 



our calculations, we first find B^ using (|6,256j ) but ignoring G n i (i.e. ignoring the nonlinear 
corrections). Substituting ( |6.256 ) into fl6,257| ), and also ignoring J- n \, we can calculate the 
free energy. At T = 0, we can alternatively use ( |6.253| ) to calculate the free energy. Numerical 
evaluation, shows that the free energy ( |6.253j ) gives rise to a centered rectangular vortex lattice, 



aligned with x or y axes, but now with the apex angle (3 < 60°; unlike (5.188) which gives a 
(3 > 60°. In order to map out the complete equilibrium H-T phase diagram we have carried 
out a numerical computation of the vortex lattice structure using the full expression for the 
response tensor Q(k), as given by Eqs. (6.243) and ( 6.252|) . We find that the free energy has 
two local minima for centered rectangular lattices aligned with two high symmetry directions 



(see Fig. 3.14), as expected from the tetragonal symmetry of the problem. Which of the two 
becomes the global minimum depends on temperature and field. The results are summarized 
in Fig. 6,15| . For high temperatures, the exact result agrees well with the one obtained from 



the long wavelength free energy (6.248). The deformation of the lattice from perfect triangular, 
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Figure 6.15: Equilibrium angle (3 as a function of reduced temperature t = T/A^ for various 
fields. Open symbols mark lattice with orientation along x or y direction while solid symbols 
mark the lattice rotated by 45°. We use A = 1400 Aand k = 68. Inset: (3 as a function of field 
at fixed T. 



grows with decreasing temperature, reaches maximum, and then falls. Maximum distortion 
occurs around t ~ 0.3, attaining (3 — 70° at 10T. Extrapolating this field dependence (see inset 
to Fig. 6.15|) , the lattice should become square lattice around H « 30T, but this field is outside 
the domain of validity of the London model. At lower temperatures the distortion decreases but 
instead of going all the way back to triangular at t*, the lattice undergoes a first order phase 
transition to another centered rectangular lattice rotated by 45° and with (3 < 60°. Further 
decrease of temperature causes the angle to grow again. We note that the precise temperature 
at which it crosses 60° depends on field, but for all fields it is close to t* = 0.19, as predicted by 
the long wavelength approximation. At yet lower temperature we predict another first order 
transition to a centered rectangular lattice along x (or y) with f3 < 60°. It has to be noted 
that the free energy difference between the two minima is very small in the region where the 
45° rotated lattice wins. It is thus likely that real system, in which vortex pinning to various 
defects occurs, will prefer to remain in the metastable state at these intermediate temperatures. 
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Figure 6.16: Apex angle [3 as a function of magnetic field B at T = 0. Circles represent 
the result of the calculation using only the nonlocal corrections. Squares correspond to the 
calculations considering both nonlinear and nonlocal corrections. 



Instead of two consecutive first order transitions the experiment would detect only a smooth 
crossover from a lattice with (3 > 60° to the one with (5 < 60°. Alternatively domains with 
various orientations may develop. 

At this point, we can add the nonlinear corrections to our numerical calculations. The 
general strategy is to use ( |6.256 ) to calculate iteratively. At each step G n i has to be 
calculated numerically. The non-analytic form of G n \, prevents us from using usual convolution 
integrals in our Fourier transformation. Instead, we first calculate G n \ on a lattice in a unit 
cell of the vortex lattice in position space. Then we use Fast Fourier Transformation (FFT) to 
calculate the Fourier transformed G n \ at the reciprocal lattice wave vectors k. In our numerical 
calculation, we use Ao = 1400 A, k = Ao/£o = 68 as before and also 7 = 2. 

At T = 0, the stable orientation of the lattice is again aligned with a and b axes. Fig. 



6.16j shows the results of our numerical calculations for f3 as a function of magnetic field. The 
upper curve (squares) is the result of combined calculation, considering both nonlocal and 
nonlinear corrections, i.e. using ( |6.256 ) and (6.257). The lower curve (circles) on the other 
hand, corresponds to taking into account only the nonlocal corrections (this curve is the same 
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Figure 6.17: The effective penetration depth as a function of the magnetic field. Circles repre- 
sent the result of the calculations considering only the nonlocal effects whereas squares are the 
result of combined calculation considering both nonlocal and nonlinear effects. 



as the t = curve in the inset of Fig. |6.15| ). As is clear from Fig. 3.16, the difference between 
the two cases is small and about one or two degrees. Therefore, the phase diagram given in 
Fig. |6.15 retains its validity qualitatively, even after adding the nonlinear correction to the 
generalized London free energy. 

We also calculate the effective penetration depth A e g- for different magnetic fields in almost 



the same way as it is calculated from /iSR data | 125 |. In these experiments, the //SR precession 
signal obtained from experiment is fit to a signal obtained by Fourier transforming a theoretical 
magnetic field distribution function n(B) defined by 

1 



n(B') = ^ / d 2 r5[B' - B(f)\ 



(6.259) 



with Q, denoting the area of a unit cell. The magnetic field B(r) in ( |6.259| ) is calculated on 
a hexagonal vortex lattice with the same average magnetic field as the experimental field and 
using the ordinary London model with some cutoff function. The A which provides the best fit 
to data is considered as A e g. Here, we calculate A e g- in a different way (but similar in spirit), 
using the fact that in the ordinary London model, for a general vortex lattice and for a large 



Chapter 6. Microscopic Model 



100 



enough field, 



(B - BY ee AB* = B 2 £ - A- 4 5 2 E ^ = r (6A- 4 (6-260) 

were r(£) is a function that depends on the structure of the vortex lattice. Associating all the 



field dependence of AB 2 in our calculation with the field dependence of an effective penetration 
depth A e ff , we can define A e fr by 



A eff 1 ABZ\ 4 

f-(xs*j (6 ' 261) 



where ABq is the mean squared value of the magnetic field (Bo(r) — B) evaluated using the 
ordinary London model on a hexagonal lattice with the same average field B and with the 
penetration depth Ao- 

Fig. |6.17 shows the result of our numerical calculation for A e g using (6.261). The lower 
curve corresponds to the calculations including only the nonlocal correction. The upper curve 
corresponds to the result of the calculations using both nonlinear and nonlocal terms. The effect 
of the nonlinear term to the field dependence of A e g is almost nothing but an overall shift. Fig. 



6.18 exhibits magnetic field distribution n(B) at the average magnetic field B = 5.9T. The 



solid line in Fig. 3.18 represents the magnetic field distribution calculated from the nonlinear- 
nonlocal London equation. The double peak feature is a sign of existing two different saddle 
points in a unit cell of the lattice which is a result of having a (3 ^ 60° . This line-shape is then 
compared with another line-shape (dashed line) obtained from an ordinary London calculation 
but with a larger value of Ao- Similar line shapes with some additional broadening can also 
be produced from /xSR data. The broadening is due to lattice disorder, interaction of muons 
with nuclear dipolar fields and finite lifetime of muons. The resolution of the magnetic field as 
a result of this broadening is 5B ~ 10~ 3 T. The difference between the solid line and dashed 



line in Fig. 6.15 as well as the double peak feature of the solid line is therefore not observable 
by /LiSR experiments because of these broadening effects. Thus as far as these line-shapes are 
concerned, it is difficult to distinguish a nonlinear-nonlocal effect from a simple shift in the 
magnetic penetration depth in the ordinary London model. Fig. 6.19| compares the effect of 
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Figure 6.18: Solid line: Magnetic field distribution obtained from the nonlinear-nonlocal 
London equation with B = 5.9T and Ao = 1400 A. Dot-dashed line: field distribution 
obtained from an ordinary London equation on a hexagonal lattice with the same B and Ao- 
Dashed line: Result of the same ordinary London calculation but with Ao = 1850 A. 

including both nonlinear and nonlocal corrections to the London equation with the effect of 
including only the nonlocal term. Comparing the two line-shapes, it is apparent that the effect 
of the nonlinear term is small compared to the nonlocal term as was emphasized before. 

6.4 Conclusion 

Ordinary London theory is not adequate to describe all the different properties of a vortex 
lattice in high T c superconductors; especially the properties resulting from the presence of 
the superconducting gap nodes or other anisotropies on the Fermi surface. However, as we 
established throughout this thesis, a generalized London model with appropriate higher order 
corrections which take into account these anisotropic effects, can still provide a fairly simple 
way to calculate different properties of a vortex lattice. The corrections we found to the London 
equation were classified into nonlocal and nonlinear ones. In both analytic and non-analytic 
(low T) cases, nonlocal corrections play the dominant role in determining the vortex lattice 
properties. 
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Figure 6.19: Solid line: Magnetic field distribution obtained from the nonlinear-nonlocal 
London equation. Dashed line: Magnetic field distribution obtained from a London equation 
including only the additional nonlocal term using the same parameters. 

Equilibrium vortex lattice geometry exhibits novel field and temperature dependence owing 
to the fourfold anisotropic effects expressed by the corrections to the London equation (Fig. 



5.11, Fig. 6.15 and Fig. 6.16|). Numerical calculation of the lattice geometry is rather insensitive 



to the details of the vortex cores. The reason is that the details of the magnetic field inside 
the vortex cores mainly affect the magnetic self-energies of the vortices. In magnetic fields 
far below H C 2 on the other hand, the vortex lattice geometry is mostly determined by the 
magnetic interaction energy between vortices which is not sensitive to the precise shape of the 
core. Therefore, our replacement of the vortex core by a simple Gaussian source term should 
be adequate for the vortex lattice structure calculations. 

The effective penetration depth A e g- also exhibits field dependence at low temperatures as 



illustrated in Fig. 6.17. Some important points need to be emphasized here: 



(i) The field dependence of X e g is not linear as is evident from Fig. |6.17 . Variation of A e g 



with magnetic field is faster at lower fields. At low magnetic fields, the relative variation of the 
effective penetration depth in our calculation is about 7% for an increase of IT in the magnetic 
field which is close to 7.3% variation obtained from /iSR data for an optimally doped and 9.5% 
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Figure 6.20: Comparison between our nonlinear-nonlocal theory and a recent /xSR result at 
T = 0. Squares represent our calculation of the effective penetration depth using nonlin- 
ear-nonlocal model (cf. Fig. 3.17| ), Circles correspond to the recent //SR data [158]. Our Ao (= 
1400 A) is a fitting parameter which might be different from the measured value. 



variation for a detwinned underdoped YBCO single crystal using the same cutoff function as 
( 5.168 ) for fitting calculations |125| , |127|] . Most recent fiSR results [158] extrapolated to T = 0, 
has provided a measurement of A e ff as a function of magnetic field up to 7.5T. An excellent 



agreement between this result and our nonlinear-nonlocal theory is presented in Fig. 6.20 



(ii) More importantly, this field dependence of A g has a predominantly nonlocal origin 
rather than a nonlinear one; contrary to what is generally believed. The contribution of the 
nonlinear term to the total (minimized) free energy is almost one order of magnitude smaller 
than the nonlocal term. What is more important however, is the field dependence and (3 
dependence of these terms not their orders of magnitude at fixed (3 and B. As we mentioned 
earlier, we consider (3 as a variational parameter to be fixed by minimizing the London free 



energy. As can be inferred from Fig. 6.16 and Fig. S.17, the field dependence and (3 dependence 



of the nonlinear term in the free energy is also smaller than the nonlocal term. It is worth noting 
that in the Meissner state, a linear nonlocal term can never produce field dependence in the 
penetration depth (as it is usually defined in that state) and therefore a nonlinear term is 
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necessary for such an effect [104]. In the vortex state on the other hand, a nonlocal term can 
result in a field dependent effective penetration depth, if we define it the way it is defined in /iSR 
experiments. We should emphasize here that in this way of definition, does not directly 
give the superffuid density as it does in the Meissner state.. 
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Figure 6.21: Magnetic field as a function of the distance from the center of the vortex 
(a-direction) for an isolated vortex. The solid line corresponds to the nonlocal London equation 
whereas the dashed line represents an ordinary London calculation with the same value of Ao- 
Bq = B(r = 0) for the ordinary London vortex. 



To understand this better, let's neglect the nonlinear term and assume a linear but nonlocal 
London equation. The total magnetic field, in this case, will be the superposition of the fields 
around individual vortices. The magnetic field around an isolated vortex is given by 

d 2 k F{k)e i% - i> 



B(f) = $ c 



(6.262) 



(2vr) 2 l + Qji^hkj 

where $o is the flux quantum and F{k) is the cutoff function resulting from the source term. 
£ij(k) is defined in ( 3.255| ). The field decays with some constant (field independent) decay rate 
A; which is the Meissner state penetration depth (see Fig. 6.21| ). Increasing the magnetic field 
does not affect the profile of the magnetic field and therefore A. Rather it squeezed the vortices 
together. For small values of k, Cij{k) Aq<%. Since small k corresponds to large r, one 
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expects an isotropic field, similar to the local London case, far away from the vortex core. For 
large values of k on the other hand, Cij(k) has strong ^-dependence with four-fold anisotropy. 
Thus the closer to the vortex core, the more deviation form an isotropic ordinary London single 



vortex is expected, as is clearly shown in Fig. S.21. At low magnetic fields, the vortices are 
far apart and their magnetic fields overlap in regions far away from their cores. The properties 
of the vortex lattice should then be similar to the ordinary london hexagonal lattice. As the 
magnetic field is increased, the vortices come closer to each other. Although the profile of the 
magnetic field around each vortex remains unchanged, the overlap regions will be closer to the 
vortex cores and will be more affected by the nonlocal term. Therefore, it is conceivable that at 
large magnetic fields, the vortex lattice properties such as the magnetic field distribution will be 
affected by the nonlocal term in the London equation. The magnetic field near the vortex core 
is always reduced by the nonlocal term as is clear in Fig. |6.21 . This is because Cijik) — Aq^- is 
a positive definite tensor for all k and therefore the denominator of the integrand in ( |6.26S ) is 



always larger than the corresponding ordinary London one. As a result, the magnitude of AB 2 
is smaller for the nonlocal case and therefore A e g- defined by ( S.261| ) tends to be larger. This 
explains why A e fj/Ao is always greater than one in Fig. |6.17 . Fig. |6.18| exhibits the resemblance 
between a change in the magnetic field distribution due to the nonlocal term and due to a shift 
in the ordinary London penetration depth. The slight difference between the solid and dashed 
lines in Fig. |6.18 would be unobservable in /j,SH experiments as a result of the broadening 
effects. Since no field dependence due to the nonlocal term is expected in the Meissner state, 
A e ff as defined here and also in /iSH experiments is expected to be conceptually different from 
what is usually defined as penetration depth in the Meissner state, although they are closely 
related. 

(iii) The magnitude and field dependence of A e g is not so sensitive to the apex angle @. In 
other words, a few degrees change in the variational parameter (3 does not modify the magnetic 
field distribution as much as a variation in the average magnetic field does. 

(iv) Calculation of A e g- is rather sensitive to the form of the vortex source term. The 
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importance of the source term in the calculation of AB 2 has already been emphasized by 
Yaouanc et al. |159| ]. In Ref.[125], A c g is obtained by fitting to the //SR data using both a 
Gaussian cutoff ( J5.168 )) and also the cutoff function proposed by Hao et a/.[16C, 158]. The 
difference between the two cases is significant and about 30% for the magnitude of A e fr and 
even more (for detwinned sample) for the relative variation with respect to the magnetic field. 
This can explain the importance of the source term in calculations of the effective penetration 
depth. 



/iSR experiments [|125| , |12g|] on NbSe2, which is believed to be a conventional superconductor, 
also show a field dependence in the effective penetration depth; although it is much weaker than 
what is observed in high T c compounds. Since there is no node in the superconducting gap 
of these materials, the theory presented in this thesis cannot explain this field dependence. 
However, since the size of the vortex core in these materials is large and comparable to the 
vortex lattice spacing for the magnetic fields of experimental interest, it is conceivable that a 



significant effect can come from the cores as is pointed out in Ref. [159]. Thus, a more careful 
consideration of the vortex core might be necessary in order to have a better quantitative 
explanation of the experimental results. 
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